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A Method and System for Discovery of Trades Between Parties 



1 Related Applications 

This application claims priority to provisional application no. .60/168,754 filed on December 6, 1999, ti- 
tled, "An E-Commerce Infrastructure for Value Chains", the contents of which are herein incorporated by 

reference. 

2 Field of the Invention 

The invention relates to a method and system for discovery of trades between parties. In particular, the 
invention is a system which allows buyers to define their preferences and sellers to define their capabilities, 
then determines which trading points maximize the utility of the buyer. The system suggests trades by 
exploiting the flexibilities and tradeoff} encoded by both parties, thus providing win-win.trades. A second 
level of optimization ranks the trades with all suppliers, allowing the buyer to rapidly determine the best 
alternatives. The system allows for rich negotiation spaces and supports continuous, discrete, and range or 
interval decision factors. 



3 Background of the Invention 

The present invention relates to methods of automatic exploration and exploitation of the flexibilities pos- 
sessed by negotiating parties to uncover improved win-win agreements. The invention describes computa- 
tionally efficient mechanisms that are applicable whether there are one or many selling parties. The precise 
number and types of negotiating dimensions are irrelevant as long as they are numerical. Thus the present 
invention applies equally to the optimal determination of terms in the purchase of a commodity or an 
arbitrarily complex artifact. 

Electronic markets have proliferated over the last few years with the advent of B2C (business-to-consumer) 
and B2B (business-to-business) electronic commerce. Such market places have yielded significant cost savings 
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by lowering the transaction costs between buyers and sellers. Buyers have also profited through increased 
^petition between suppliers. However, electronic markets have hurt suppliers, since the zero-sum neg- 
ation over price has been at their expense. The present invention describes a tool whereby cost savings for 
both parties are derived from the discovery of win-win trades. Mentally, the system works by allowing 
trading parties to describe their desired trade across multiple dimensions and to express their flexibihty 
around this ideal trade. Through an algorithmic exploration of their flexibilities, the present invent™ can 
discover trades that are near the ideal trades of both parties, enabling both to win. • 

The adoption of B2B and B2C electronic commerce was facilitated by the migration of catalogues onhne. 
This familiar method of presentation ameliorated the significant cultural change to electronic trade. For the 
foreseeable future, electronic commerce will be dominated by online catalogs. At present, online catalogues 
are direct translations of their hardcopy counterparts where the attributes of a product are described and a 
price quoted. Inevitably however, online catalogs will become more expressive. Catalog entries will be able 
to represent price breaks for large quantity orders, lot sizes, etc. Thus it is important that any software 
(like the present invention) that uncovers mutually beneficial trading scenarios is able to operate with such 
catalogs. Consequently, in the present invention there is an asymmetry between buyer (usually a human) 

and seller (usually an online catalog). 

One of the reasons catalogs have come to dominate electronic commerce is that the types of goods that 
can be represented in catalogs are simple. Whether the product is pens or paper dips, different vendor's 
offers differ little from each other (a pen is a pen is a pen), and a quick scan of a catalog gives a buyer 
enough information to make an informed purchase: These types of goods are low margin and inexpensive. In 
contrast, the vast amount of purchasing between businesses involves materials which are directly connected 
with business operations - car parts, turbines, etc. Such direct goods are the future of electronic commerce 
Unlike presmt-day engines, any truly useful procurement tool must be able to support direct materials w,th 
complex attributes and complex inter-relationships between its components. 

Electronic commerce offers unprecedented opportunities for more informed decision-making for both 
buyers and sellers. The past few decades have seen the widespread adoption of enterprise resource planning 
(ERP) systems, to the point that now almost every major company has some form of ERP software. ERP 
functions as the digital nervous system of a company, transmitting and logging information between the 
company's many different business functions. ERP software keeps track of inventory, monitors the state 
of purchase orders, signals when a company should reorder direct and indirect materials, and a mynad of 
other functions. Consequently, EBP databases are a rich source of information to optimize a company s 
operations. Yet today this information is rarely used t» make more informed buying and selling demons. 
The present inventioncan utilize such formation sources to optimize a company's interactions with supphers 

and customers. 
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One important manner in which this optimization can occur is through an analysis of oil cost factors. 
Current buying and selling practices often focus on limited goals, e.g., minimize the total purchase price. 
Myopic purchasing strategies often result in higher total cost of ownership when all cost factors relevant to a 
product in its lifetime of use are included. These other cost factors can be significant. Why save the money 
in taking delivery two days late if the receiving docks will be full at that time and an additional shift needs 
to be hired to clear the docks? Why order the cheaper drill bit if it is much more expensive to replace when 
it breaks? The present invention improves trades by minimizing the total cost of ownership of a product, 
yieldingsignificantsavingsto its users. Many total cost factors are difficult to quantify - e.g. what is the cost 
of dealing with a unionized versus a non-unionized supplier? Consequently, the present invention supports 
qualitative (best guesses and intuition) as well quantitative factors. 

All companies are situated in a supply or value chain. At each step in the chain, a company purchases 
from its suppliers, transforms these inputs, and sells the output to its customers. The termination of the 
supply chain is the sale of the final product to the end consumer. Since the only influx of external capital 
comes from the end consumer, companies have realized that they compete not only as individuals but 
also as entire supply chains. As result, software products have recently become available which attempt 
to streamline the operations of links within the entire supply chain. This software, variously called supply 
chain optimization (SCO) or advanced planning and optimization (APO), operates on the basis of forecasted 
demand at various points within the supply chain. Based on these predictions, plans are generated telling 
companies how much to produce and how to schedule their operations. SCO systems are a valuable source 
of intra^ompany information - data the present invention capitalizes on. Because SCO software relies on 
forecasted demand, it is only as helpful as the forecast is accurate, and, unfortunately, in many cases demand 
is very difficult to predict. How can the software know that laundry detergent will go on special at grocery 
stores in the Northeast in 7 weeks? As a result of the volatility in demand and the many other unpredictable 
perturbations that plague supply chains, companies keep significant buffers in the form of inventories. In 
addition to planning, businesses must also be able to adapt to unplanned effects. Such adaptation requires 
flexibility and a means to exploit that flexibility. The present invention exploits the flexibility of trading 
parties to streamline the operations of supply chains by smoothing the boundaries between trading parties. 

The present invention is therefore a system to allow trading parties to express trading desires and con- 
straints across many and varied different factors. These trading preferences are informed by many different 
data sources to optimize for a company's internal operations and its connections to it's supply chain through 
an analysis including total cost factors. The- flexibility expressed by all trading parties is exploited to locate 
win-win opportunities for all parties if they exist. 
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4 Summary of the Invention 

We describe the present invention in its application to facilitating trade between buyers and sellers, but 
note that the mechanisms described are much more general. We can easily imagine, for example, using the 
present invention to match individuals (with the desires and skills) to projects. 

The inspiration for the present invention comes from utility theory developed by economists since the 
1960's. Since we are interested in multiple dimensions of negotiation, we draw from the multi-attribute 
utility theory literature. 1 Utility is an abstract concept which has been formalized in various ways. For the 
present purposes utility, u, is a number between 0 and 1 representing a party's willingness to trade. Larger 
values indicate a greater willingness. 

4.1 The Negotiation Space 

In any negotiation the parties must come to agreement on the factors requiring negotiation. We call these 
factors dimensions or variables. As an example, when purchasing a car, the buyer may be concerned with 
price, time of delivery, and color. Each factor price, time, and color is a dimension. Most dimensions can 
be classed as one of three types: continuous, discrete, or range/interval. A continuous dimension is one like 
price for which the buyer's utility varies smoothly across that dimension. The buyer's utility at $23 001.00 is 
almost the same as the utility at $23 000. Color is a discrete dimension. Since the car may only be available 
in black, white, and silver, the domain of this dimension is the finite set of values {black, white, silver}. 
Moreover, the buyer's utility may be quite different for the three colors. The third class of dimensions is 
called interval dimensions. An interval dimension arises often in B2B negotiations. H a machined part is 
built to some tolerance (e.g., the inner diameter of a screw is between 24.5 and 25.5 mm), the range of 
variability in the dimension is specified as an interval. In the language of statistical quality control, a certain 
percentage of the machined parts will fall in this range. These three broad classes of variables capture almost 
all the types of attributes relevant to B2B negotiation. 

The present invention operates over any number of continuous, discrete, and range or interval variables. 
We call the negotiation space X and any point in the negotiation space (x,x,r) 6 X. It is important to 
recognize that the single trading point (x,x,r) may have multiple components, e.g., price = $23 000, time 
of delivery = 3 weeks, color = black. 

In the present invention, the space of negotiation is agreed upon by all parties involved prior to the 
commencement of any negotiation. We can, however, imagine more dynamic situations in which dimensions 

are introduced and discarded o vertime. 

•For a good introduction to multi-attribute utility theory see [1]. 
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4.2 The Buyer's Utility Function 

A party defines it's utility function over this space so that every (x,x,r) is assigned a utility number 
indicating the party's willingness to trade. We indicate the utility function as «((x,x,r)). A great deal 
of work has been done on the appropriate form for utility functions. In the present invention, we take 
a simple form for the utility function for two reasons. First, we would like the form of the utility to 
be conducive to rapid computation. Second we would like the utility to be simple enough to be easily 
understood by and elicited from users of the invention. With no loss in generality, we write the utility 
function as u((x,x,r)) = exp(-d((x,x,r))) where d(x) is interpreted as the distance of trading point 
(x, x, r) from the most preferred trade. 

So that we can operate against seller catalogs, only the buyer needs to define a utility function. Across the 
continuous dimensions, the buyer's utility is denned by specifying the most preferred (or ideal) continuous 
dimensions and the manner in which utility drops off as we move away from this ideal. For the discrete 
dimensions, the utility is specified in tabular form since there are a finite number of alternatives. Again, the 
buyer must specify it's ideal discrete values and how utility decays away from those values. In section 6.1 we 
describe how this is accomplished. The range dimensions contribute to utility similarly; the buyer specifies 
an ideal range and the utility decays for ranges other than the ideal according to their distance from the 
ideal. 

The utility function can also express tradeoffs between variables, e.g., I may take delivery in 5 weeks if 
the price drops to $20 000, or I may accept the white car if I can take delivery in 2 weeks. The tradeoffs 
may be between pairs of continuous dimensions (as in the first case), between pairs of discrete variables, or 
between continuous and discrete variables (as in the second case). 

4.2.1 Normalization and Weighting 

When utility is defined over different types of variables, it is important to normalize the contributions of 
each variable so that the buyer can weight the importance of the various contributions to utility. This is 
a difficult problem. How should a buyer's color preferences be normalized so that they can be traded off 
against time of delivery? The present invention solves this problem by requiring that the average distance of 
any negotiation variable from its ideal value is the same for all dimensions. Since the buyer is more interested 
in those regions of the negotiation space where the utility is high, the average is weighted by utility. This 
procedure defines a manner in which to define a baseline where all dimensions contribute equally. Given this 
baseline, the buyer can then weight the various contributions and obtain useful results. 
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4,2.2 Utility Elicit ation 

Since utility is fundamental to the present invention, its elidtation from the buyer is important. Utility may 
be defined using any of a number of sources: 

1. graphical user interfaces associated with the invention 

2. standard benchmark criteria applicable to the domain 

3. formal methodologies like the analytical hierarchical process [2], or discrete choice analysis [3] 

4. inferred through models 

We expand briefly upon method 4. As discussed in the background section, it is important to buyers to 
minimize their total cost of ownership. If we have a function representing these costs as a function of the 
negotiation variables, and perhaps other factors, this function can be used to infer a utility function which 
will act to minimize the total costs. Later we describe how this can be accomplished. 

4.3 A Supplier's Capabilities 

As noted previously, suppliers are treated differently from buyers so that the tool can operate against catalog 
information with no human intervention required on the part of the seller. In fact, we do not require sellers 
to define a utility at all. 

A supplier cannot offer deals at all points within the negotiation space X, e.g., he certainly can't offer 
the black car tomorrow for free! A capability then represents the ability of a supplier to deliver and defines 
a subspace of X. It can include such things as price discounts on large volume orders, variation in delivery 
time as a function of price, etc. Since these relationships are already specified by businesses in terms of 
simple rules like "the price per unit i$ $10.00 if 1 to 999 units are ordered and $9.50 per unit if 1000 or 
more units are ordered n } suppliers' capabilities are represented in the present invention by piecewise linear 
functions. 

4.4 Negotiation Constraints 

Both parties may have constraints which must be satisfied in order for them to trade. For example, the buyer 
may not buy the car unless he gets it within 6 weeks, or he may not purchase the car if it is available only in 
white. These are examples of continuous and discrete constraints, respectively. A continuous constraint sets 
a requirement on the continuous variables. In the present invention, continuous constraints must be either 
linear or quadratic. Discrete constraints involve discrete variables. A discrete constraint can be expressed 
as a list of the allowed (or disallowed) combinations of the discrete variables for which the trade will be 
acceptable. For example, if the buyer would accept either the black or the silver car, the constraint would 
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B* both these colors as viable. It is important to note that both continuous and discrete constraints may 
involve one or more variable.. We can also express constraints involving both types of variables by allowmg 
the continuous constraints to differ depending on the discrete variables. 

4.5 Utility Optimization 

With the major components of the invention in place, we describe how the overall system works. As a 
procurement tool for the buyer, there are two levels of optimization. First, for any given supplier we 
maximize the buyer's utility, subject to the supplier's capabilities to find that trade which makes the buyer 
as happy as possible. Since we are optimizing within a supplier's capabilities, the supplier has expressed a 
willingness to complete the trade at whatever point is determined to be optimal. Hie tool then optimizes 
across suppliers to rank them according to utility at the optimal point. A graphical user interface allows a 
buyer to investigate the trades suggested by the tool by sorting according to any dimension or by the overall 

utility. • % a 

Utility, while a useful concept in assessing an overall score, may be of limited use to a buyer due to * 
abstract meaning. Consequently, we can also applythe total cost of ownership function to the result to rank 
order the suggested trades according to their various cost components. Recall that for any trade x 6 X, the 
total cost of ownership function returns the various cost contributions. This additional information aids the 
buyer in his purchasing decision. The utility number for each trade is still useful because the total cost of 
purchase function includes only those cost factors which can be quantified, whereas the utility also mcludes 
"softer" qualitative factors. 

4.5.1 Aggregation 

In addition to optimizing against one supplier at a time, the present invention can also be used to optimize 
against an arbitrary aggregation of suppliers. This is important if, for example, no single seller can supply 
the large volume requested by a buyer. In this mode of operation, the buyer specifies sets of supphers 
participating in the aggregation and the dimensions over which aggregation can occur, and the tool finds 
the optimal combination in which to distribute the volume dimension over the allowed suppliers. 

5 Brief Description of the Figures 

Figure 1 shows an an architecture for the invention. 

Figure 2 shows a schematic of a buyer-specific capability with examples indicating potential input. 
Figure 3 shows a schematic of a suppUer-specific preference with examples indicating potential input. 
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n c 
rid 



number of continuous dimensions 
number of discrete dimensions 
number of range dimensions 



n c - vector of values for continuous dimensions 
nd-vector of values for continuous dimensions 



value of ith continuous variable 
value of ith discrete variable 



Table 1: Definition of the negotiation search space. 

6 Detailed Description 
6.1 Theory 

In this section we outline the mathematical foundations of the optimization process in sufficient detail to 
allow for computer implementation. 

6.1.1 The Negotiation Space 

In Table 1 we define the parameters which collectively define the space of negotiation X . For each of the 
n c continuous variables, we specify an allowed range over which that continuous dimension may vary as 
Xi£Xi = where x is the n c -vector of lower continuous bounds and * is the n c -vector of upper 

continuous bounds. Each discrete variable assumes a value from within its domain x { € 2>*. Without loss 
of generality, we label the domain of discrete variable i by V { = [1, ■ • ■ . *] where 4 > 0 is an integer giving 
the number of possible values discrete variable Xi may assume. 

With these definitions, we define the space of negotiation by the tensor product X = X x 0 • • • ® X« c ® 
Vi ® • • • ® X>„ a . Range variables are treated separately and not negotiated over. 

6.1.2 The Utility Function 

The utility function is a mapping from X into the interval [0, 1]. As indicated earlier we assume the utility 
to have the form u(x,*) = exp[-d(x,x)] where d(x,x) is interpreted as a distance. In what follows we will 
assume that in its simplest form the distance function has the form 

d(x,x,r) = (x-Mx^Mtx-^)) +Z(x) + J R(r;x). 

Each contribution to the distance function is positive. We consider each contribution to the distance in turn, 
beginning with the range variable contribution ii(r;x). 
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First, we note that the range distance depends on the setting of the discrete variables. This allows 
the buyer to express different preferences for the range variables depending on discrete factors. The total 
range distance is summed up over all possible range variables so that Ji(r; x) = ££i Jfc(r i; x). The vector r 
indicates the preferred values for all range variables. If range variable i is specified as the interval n = (u>?i) 
(where r< > ij) then r is an n r -vector of such tuples. The distance contribution, IU, from one range variable 
will depend on the application. If the range variables are meant to represent the tolerances on machined 
parts where issues of statistical quality control are important, then the distance between two ranges might 
be related to the overlap between Gaussian distributions. If n is interpreted as a Gaussian having mean 

+ ?i)/2 and standard deviation proportional to n - n { then an appropriate range distance is given in 
Appendix A. Other choices for the range distance function are certainly possible. 

The continuous distance is quadratic and determined by the positive semidefinite n c x n e matrix C _1 . 
We have allowed this matrix to vary with the setting of the discrete variables and indicated this explicitly 
through C -1 (*)- a The «c-vector /i may also depend on x and indicates the point at which the utility is 
maximal - p is thus identified with the ideal value for the continuous variables. The precise quadratic form 
is convenient, but, using recent developments in interior point methods, other convex functions are also 
computationally tractable [4]. 

The discrete distance is determined through the function Z[x) which maps the discrete space X>i ® • • • 0 
V ni onto the positive real line [0,oo]. In keeping with the assumption that distance is a function of only 
pairs of components x h Xj , we assume the discrete distance has the form 3 



i=i i=i(^o 

Each contribution Z id is a table consisting of (fcrf, entries, where ^(xi.xj) can be interpreted as the 
distance if discrete dimension i has value xj conditioned on discrete dimension j having value xj. The 
diagonal terms Z< offer an unconditional distance. The most preferred value for the ith discrete dimension 

is that for which Z*(x*) = 0. 

Rather than require the user to enter the distances explicitly, there are numerous ways in which the 
distances can be generated automatically based upon a buyer's ranking of preferred values. Further details 
can be found in Appendix B. 
Weighting of Dimensions 

In many cases it is important for simple modifications of the distance function to re-weight the contributions 
to the total distance. If w c is an n c -vector of weights for the continuous dimensions, we can accomplish this by 
letting C" 1 = W e C- l W c where W c is the diagonal matrix W e = diag(w c ). 4 In a similar way we modify the 

2Where no confusion will arise we eliminate the dependence of C _1 on x for notational simplicity. 

3 Later -we shall generalize this distance to include weighting of dimensions. 

4 M = {mi j) = diag(v) is the diagonal matrix formed by setting m.i = »i and "WJ = 0 for 3 £ »"• 
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discrete distance to Z^jfr, *,) = «n,«w«Su(**. *i) where w <* is **» «<i-vector of weights for the discrete 
variables and w^t is its tth component. The range contribution is also modified so that i^fa) = wtsiJfafo) 
where w r is the n r -vector of weights for the range variables and w rii is its ith component. For convenience the 
weights are normalized so that (l«w e ) 2 +(l'w d ) 8 +l'w P = 1. With little additional complexity the dimension 
weights can be made dependent on the setting of the discrete variables but we will assume throughout that 
the weights are constant. 

With these modifications, the total distance function becomes 

d(x,x) = (x-rtx^C^xXx-Mx)) + Z„(x) + iUr). . 0) 

where Z*{x) = E2i «*«{w««Z«(*») + and J^(r) = E£i 

Assigning weighting factors is useful only if the relevant contributions have been previously normalized 

so that they are all roughly the same magnitude. This serves as the baseline for which all weights are equal. 

The question immediately arises as to what criteria to use to weight the distance contributions. 

We shall determine scaling factors, Q d > 0 and Q r > 0, so that the average distances per dimension of 

the discrete, range, and continuous contributions are equal, where by average we mean a utility-weighted 

average over the entire space of possible trades. This weighting places more emphasis on the better trades 
If d c , dd, and dr are * te continuous, discrete, and range contributions to the total distance, then after 

multiplication by the scaling factors d = d c + Q d d d + Q r dr. The scaling factors are determined through the 

utility weighted average distances defined by 

. . _ E~E,/vrf"Qrdrexp[-Q r d r -Qddd-d e ] g M exp[-Q, 1 4,]E r d r exp[-3 r d r ] J v du exp[-d e ] 
(dr> = E~Er Sv du ecPl-Grdr - Qdi - del Wr Ex «rf-Qd*<| Er exp[-Qrdr] Sv du exp[-d c ] 

_ ^^LrSydviQitUexpi-Qrdr-Qddi-dc) ^^[-QM E r exrf-Qr^l Sv du exp[-d c ) 
{dd) = Ex Er Sv du expl-Qrd, - Q d d d - d c J Vd Ex e*v\-Qdd*] Er exp[-<?rdr] f v du exp[-d c ] 
. ExErJv<fa<* e exp[-Q r dr- Q d d d -d c ) Exexp[-Qdd j ]E r exp[-Q r d r ]/ v dud c exp[-d c ] 

{de> = Ev Er Sv dU «Pl"Ordr " Qddd - rfcj Ex «p[-Qd*J Er exp[-Q r dr] S V da exp[-d c ] 



A few comments on the above equations are in order. First, Ex indicates the repeated sum E 
over all possible discrete trades. £ r indicates a sum over all the range variables and the integral over 
volume V indicates integration over the continuous trading volume of interest. Finally, we have not included 
a scaling factor Q c on the continuous distance, since this can be made equal to 1 if we reinterpret Q r as 
Qr/Qc and Qd as Qd/Qc* Each of the averages is an explicit function Q d and Q r . 

The requirement on equal average contributions determines the two unknowns Q r and Qd through the 
equations: (rf r )/n r = (dc)/n c and {d d )/n d = (d c )/n c . These two nonlinear equations are coupled in terms 
of Q T and Q d and must be solved simultaneously for Q T and Q d . Rirther details are found in Appendix C. 

6 We 8 i Q gted the continuous variables out, since there will always be continuous variables in any trading scenario. 
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Table 2: An example set of constraints involving 3 variables where the domains of all variables are V\ = 
X> 2 = X> 3 = {1,2,3}. Constraint (a) requires that the values assumed by xi, x 2 , and x 3 are all different 
from each other, and constraint (b) requires that the value assumed by xa is even. See text for the solution 
to both constraints. 

6.1.3 Constraint Specification 

Buyers and sellers may express constraints over both continuous and discrete variables. 
Continuous Constraints 

For simplicity (and because additional expressiveness is rarely required) we assume that the buyer's con- 
straints over the continuous variables are linear. 6 This allows a buyer to express a constraint, e.g., the time 
of delivery must be within 10 days or I will not trade, i.e., t < 10. We allow for both inequality and equality 
constraints which can be expressed as G^x = and G^'x < If there are dp equality constraints 
then has c\ } rows. Similarly, has 4 rows if there are q 1 inequality constraints. We allow the 
constraints to depend on the setting of the discrete variables, and to be explicit we often write G^(x), 
g< 6) (x), G«(x), andgW(x). 

Discrete Constraints 

We use a standard methodology to represent and process constraints over discrete variables [5]. Abstractly, 
a constraint over a (perhaps proper) 7 subset of the discrete variables is represented as a list of all the allowed 
combinations the variables may assume. An example representation of a pair of discrete constraints is given 
in Tible 2. There are two solutions to this set of constraints: x\ — 1, *a = 2, h% = 3 and x\ = 3, 
K2 = 2, xfc = 1. We indicate these solutions as {(*i,l),(*fc>2)},(x3,3)}} and {(*i,3),(xi,2)},(x3»l)}}. 
respectively. Each solution where all the variables have been identified with specific values from their domains 
is called a labelling. 

6 With little additional complexity we can also handle quadratic constraints, see [4]. 
7 A proper subset of a set is the set itself. 
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Computationally efficient representations are used to ensure that only feasible combinations of variables 
processed. Numerous third-party libraries offer constraint programming functionality. 8 



are ever 



6,1.4 Utility and Total Cost of Ownership 

The buyer's utility function and associated constraints may be difficult for many users to define. In this 
section we show how models of the buyer's business can be used to define utility in a natural manner. 

We imagine a function which provides an estimate of the total cost of ownership for any given purchase. 
Cost contributions to this function might include piece part costs, freight costs, setup costs, quality assurance 
costs, repair costs, etc. It is important to include all quantifiable costs associated with the lifetime of use of 
the purchased product because it is this function we will be minimizing. Significant savings may be obtained 
by taking a longer-term view of the purchase. Revenues (negative costs) generated from the purchase are 
also included in the function so that the function represents some measure of profitability associated with 
the purchase. We write the total cost of ownership function as C o (x,x,r;0). We explicitly indicate the 
dependence on the negotiated trade parameters x, x, and r, as well as other factors 0. The other factors 
might include forecasted demand, current inventory levels, etc. These factors will vary over time, and they 
can be extracted from the buyer's ERP and supply chain management systems (SCM) in real-time just 
before the purchase to ensure continuous real-time optimization. See section 6.2.1 for further details. 

Minimization of C o (x,x,r;0) defines an ideal trade dependent on current conditions: Xo pt (0), *o P t(0)> 
r op t(0)- If desired, these can be used to define /i = Xo pt 0?, r = r opt (£) and the desired ideal discrete 
configuration x opt (0) (having distance contribution Z = 0). Moreover, the tradeoffs between continuous 
dimensions around this minimum can be obtained through calculation of the Hessian matrix H — [hij] 
where the i, j matrix element is given by 



dxidxj 



=Xo P t(0),*»<c.pt(/&),r=ropt(0) . 

We then identify CT 1 with H. In this way, little trading flexibility is obtained in directions where total cost 
of ownership rises rapidly, while significant flexibility is obtained in directions where total cost of ownership 
increases slowly. 

In summary, a total cost of ownership model defines both the most preferred trade parameters and the 
flexibility possessed around the preferred trade. The model pulls dynamically from real-time data sources to 
provide the most up-to-date optimization based on total costs of ownership and other important qualitative 
factors the buyer may wish to describe in the utility function. The same function and its constituent costs 
may also be used to help analyze proposed trades from suppliers. 

8 Sce for example vw.mozart-oz.org or wv.ilog.com. 
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6.1.5 Supplier Capabilities 

As discussed in the introduction, suppliers represent their capabilities through a specification of the subspace 
of X in which they will trade. A supplier's capabilities must specify the allowed continuous, discrete, and 
range variables. The allowed range variables are expressed as the pairs (r^r,), one for each range variable. 
For example, if a supplier produces 25mm inner diameter screws to within a tolerance of 0.5 mm, then the 
range variable is simply (24.5,25.5). These are compared with the buyer's ideal range and contribute to the 
distance function through the R«,(*) function. 

Capabilities over continuous and discrete variables are more complex. 

Continuous Capabilities 

Continuous capabilities are viewed naturally as responses to a buyer's request. Thus we distinguish between a 
buyer's requested continuous vector x< 6 > and a seller's response x<*> . A vector-valued function, f (x<« , x W , x) 
returns the response based on the buyer's request and also, perhaps, other previously defined supplier 
responses. Component /< of f defines the tth continuous variable, i.e. x\ a) = /i(x( 6 ',xW). 

Currently, suppliers are used to quoting price discounts for large volume orders and these price discounts 
are expressed as piecewise linear functions. Consequently, we restrict /* to have the following form (where 
we distinguish between the functions depending on the buyer and seller variables): 

An example of how this may be used to define a supplier response is the following: We assume three 
continuous dimensions - price, volume, and time of delivery and indicate these as [x x , x 2 , ar 3 ] = b> v » *]■ Then 
a response may be formed as 

tW=/ t (t,Wt«) = /,><■>) 

pW = /,(««,«(•>) = /„><•>) + / P(t (iW). 

The /„,» function returns the volume a supplier will fulfill as a function of what the buyer asked for. If the 
supplier can deliver any volume, this will be the identity function. K the supplier delivers only in certain lot 
sizes, this function may have a staircase shape, etc. The f t , v function indicates the time it will take a supplier 
to deliver a certain volume. So, for example, if larger shipments require longer transportation, then this 
dependence is given by this function. Finally, we turn to the price determination. In this example the price 
depends on the quantity v« being shipped and the f PtV might represent price discounts for large volume 
orders. There is also an incremental price contribution based on the time of delivery. If faster delivery is 
more expensive, this is reflected in f Pl t. 
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For a given setting of the discrete variables, each /S ) (4' ) .* W ) 811(1 /St*?*.**'*)* a one-dimensional 
piecewise linear function. Consequently, the functions can be specified by listing the breakpoints. If 
/<;>(4»>,:*W) has *g breakpoints, then we list these as {a^J = = 2 >."' = 

,$)}• and function values at these breakpoints are{^(xW(6W = l)),/g(4 B) (i# = 2 »'- ./S^SW 
k\'1))}. Similarly, is a piecewise linear function defined by the k$ breakpoints {xj^ (6$ = 1 )>** ) ( i i2 = 
2),... ,«? J (6g = and function values at these breakpoints {/g^flS = ifl./ffiO^S = 

2 )) » • • • » /i * (4* } (*>u = 4,1!)) ) • The break P° mts m indexed by the integers b$ and 6$ . 

An interval is specified by assigning a value b\'l € [1, - 1] and 6$ € [1, k$ - *] «> t^ 10 

4 ,) (biS)<4' ) <4' ) (^+D v< and 4 6) (0^4 6) <4 6) (^S+i) v-t 0) 

Within each interval the functions are linear, so we have 

4* = «(*<•>, {*B>. {©) + E-S^- W + "®>S)4 6) 

where Ci(v<'\ {6$}) « £* c| J (*<•>, ftjj )+cg(6jJ). In the above equations, the intercepts and slopes 
are given explicitly by 



and 



(.wo, _ + D/S(«P(M - «P(*a)i& ) (4 ,| (ftM 
Ci ' fc(b| ' fc) " ' 4-)(6|;i + i)-^(6, fc) ~ 



WftW , _ ^ ) (4' ) (^+i))-^ ) (4°(k0) 



respectively. An analogous result holds for the c$(b$) and m$(6$). 

To eliminate any cyclic dependence on x\ 9 ^ we must impose an ordering on x\ 8 ^ so that x\** can only 
depend on where j < L Consequently, we can write 

4'> {$.■•• MU^n^E^i^W 

Written as a matrix equation, the above becomes 

(I - M«(*M {6$}))x<»> = «(*«, {&«}, {6g}) + MW({6g})xW 

wherec(x('),{6g},{tg}) = [^.flB hflS}) - ^.{^^{bg})] 4 ,^) = [z« *«) 

°The breakpoints are assumed to be in increasing order. 

10 Note that many choices for {6$} or {b\ b J} would be inconsistent with these constraints. A simple way to fix this is to 
define a union of all breakpoints (s) and (6), and have all fa have the same breakpoints. 
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(») 



M<*)(««'),{*i2}) = 



0 
0 



n4J(xM,lB) 



0 
0 
0 



and 



M«({&g}) = 



«8<«) -MS) 
"4Mb 



In most cases x« will depend only on a subset of the variables in x<». If xW depends on n' < n of the xW 
variables, then M< 6 > is an n x n' matrix. In the example given everything depending only upon the volume 
the buyer requested. 

Since MW(xW) is lower triangular and can be inverted in time 0(n), we can rapidly express xW as 

as long as the fcg are chosen to also satisfy ^(ftg) < xj? <*1' } (*8 + Thcse C 008 ^ 0 * 8 ^ be used 
in section 6.1.6 which formulates the optimization problem. 

We also allow a supplier to express additional linear constraints so that, for example, he may represent 
that he does not deliver on Sunday. Thus the supplier may define the matrices G[ b) (x), g[ $) (m), and the 
vectors gM(x), such that G«x« = g« and G^x < g<»>. G["(x) and G<<>(*) have c« and 

rows respectively. 
Discrete Capabilities 

It is easy to imagine that a supplier's response on a discrete dimension is highly constrained by the values 
of the response on other dimensions, e.g., certain product characteristics come only in certain colors and 
package sizes. Consequently, it is not suitable to explicitly define a response but only to make available 
a supplier's constraints amongst the discrete variables. Consider 3 discrete dimensions where xi € V\ = 
KM. x 2 € 1>2 = [AB,C7,I>], and x 3 € Z> 3 - = and assume the supplier has the following 3 

constraints 

Ci(x!,x 3 ) = {{*>*)MP),{bMc>P)h 
C 2 {x* 7 >*) = {{AMB,y) 1 {D J f!)}, 

Cs(xi) = {6,c}. 
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We first note that there are 4 feasible solutions (or product configurations the supplier can meet) : [*i , x 2 , x 3 ] - 
[M»0]> [b>Di0\, [c,A,0]> or [c,A0J. Feasible solutions to the constraints define the response x<*> for the 
discrete variables. 

We indicate a supplier's or buyer's collective set of discrete constraints by C (a) (x) and C< fr > (x) respectively. 
6.1.6 The Optimization Problem 

Having defined the necessary components, we now define the optimization task which determines the con- 
tinuous x* and discrete x* parameters of the trade. 

Since the trade must be acceptable to the supplier, we maximize the buyer's utility over a supplier's 
capabilities. Equivalently, we minimize the distance from the buyer's ideal values as 

[x',x*l = arg min {(x« - M^)))*^ 1 ^)^ -/i(x«)) +Q d Z w {x&) +Q r R w W} 

where 

x« = (I - MW(x w ))" l c(xW) + (I - mWCxWJ^mWxW 

subject to the constraints over continuous variables 

d (x«)x<'> = gl (x<*>), G 2 (xM)x« < g 2 (x«) 

and the constraints over the discrete variables C< fr >(v«), CW(v«). In the above, we have defined the 
(c^ -b c^) x n c and (4° + 4 6) ) * "c matrices G^x*')) and G 2 (x«) by 



Gi(xW) = 



G^xW) 
G< 6) (xW) 



, and G 2 (x< ff >) = 



G^x^) 



The + c[ 6) )- and (c^ + c^ } )-vectors gi(xW) and g 2 (xW) are defined by 



gi 



(*») = 



g«(xW) 
gi 6) (^>) 



, and g 2 (x w ) = 



gM(*M) 



The optimization is accomplished by iterating two distinct phases. Phase one sets the continuous param- 
eters optimally for a given setting of the discrete variables. We define the functions 

d 1 (x,x) = (x-p(x))*C- 1 (x))(x-/i(x)) + fl(r;x) and cfc(x) = QdZ w (x«), 

The first phase of the optimization is the continuous problem: 11 

x*(x) = argmin di(xW,x) subjectto Gi(x)x = gl (x) and G 2 (x)x < g 2 (x). (5) 

x<-> 



u No optimization Is required over range variables, since these are specified up front by both buyer and seller and merely add 
to the total distance. 
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A detailed discussion on the solution of the phase 1 optimization problem is found in appendix D. The 
second phase determines the best value of the discrete variables by minimizing over a function of x alone 

tc* = arg min di(x(x),x) + <fe(x) subject to C< 6 >(x)AC<'>(x). ( 6 ) 

Further details on the phase 2 optimization are given in Appendix E. Once x* has been determined, we find 
x* as x* =s x(x*). 

6.1.7 Aggregation 

Often a buyer may be willing to divide an order between multiple suppliers in order to aggregate the required 
demand or to obtain better deals. In this section, we detail how the present invention supports this aggregate 
optimization. 

Aggregation can only occur over the continuous variables where values may be subdivided. Each contin- 
uous variable %i must be parcelled out amongst a set of suppliers. Consequently, we extend our notation to 
Xi -> Zi, k giving the contribution of the Jfeth supplier to continuous dimension i. The *th supplier may come 
from a (perhaps proper) subset of all suppliers. We indicate the set of potentially contributing suppliers as 
K and the number of potentially contributing suppliers as |/C|. 12 

We restrict the discrete variables to be the same across all potentially aggregated suppliers, Le., we do 
not generalize x* -> xs.fc- This simplifying assumption is made for two reasons. First, the size of the discrete 
optimization problem is smaller and so optimization be performed faster. Second, it may be difficult to 
elicit from the buyer the allowed discrete alternatives for each supplier. Nevertheless, this generalization is 
straightforward should the need arise. 

This simplifying assumption requires that the union of discrete supplier constraints C/c(x) = A*eJC C fc #) ) 
yields a feasible solution when combined with the buyer's discrete constraints CW(*)» A necessary (but not 
sufficient) condition for satisfaction is then that each constraint satisfaction problem k having constraints 
C W (x) AC£ ff) (x) has a feasible solution. 13 Henceforth, we will assume that the set of suppliers, K 7 satisfies this 
condition. If not, those suppliers violating the constraints C<*> (x) ac£ } (x) are eliminated from consideration 
in K. 

Discrete Search 

We must search over the subsets of K for feasible solutions, which is a combinatorial problem. Fortunately, 
given a complete labelling of variables, determining the largest subset is easy. For any given labelling of all 
discrete variables, if each C< 6) AC£° V k G k C K is satisfiable, then the union C w ACi** where d° = A* e * 4 #) 
is also satisfiable under the same labelling. The largest subset of variables is found by adding all k which 
13 All work to this point is thus seen as the special case \K\ = 1. 

13 It is easily seen this requirement is not sufficient by having C <*> = , 1), (xi f 2)}, c[ 9) = {(hi , 1)}, and C/ = {<*i , 2)}. 
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have feasible solutions with the buyer. We needn't worry about smaller subsets because the continuous 
optimization will assign zero values to those if appropriate. Consequently, for any given labelling x we let 
/C(x) represent the maximal subset of suppliers for which C< 6 >(x) A C/c(x) is satisfiable. It is this set of 
suppliers which enter into the continuous optimization. The number of participating suppliers is denoted by 
|/C(x)|. 

Continuous Optimization 

Optimization over the continuous variables is carried for each labelling x. Generally speaking, the buyer's 
utility will not be an explicit function of but only of a*. We assume a linear relationship between these 
two quantities so that 14 

x = Ex. 

The n c |/C(x)| vector x is defined as x* = [xi,-- ,x n J where xj = ,**|ic(*)|l- The n c x n c |/C(x)| 

matrix S is assumed known and typically has the form 16 

'i\ o* o*" 

• • , • 

where 0 is the If -vector of all zeros and & is the linear combination relating x\ to the Under our 
assumptions for E, x t - = {jx*. Li cases where the buyer wants to accumulate the results from suppliers (e.g., 
aggregating quantities) £ = 1 is the |/C(x)|-vector of all Is. In other cases the buyer may take f = 1/|/C(x)| 
so that the time of delivery becomes the average time of delivery across the suppliers. Constraints on x 
become constraints on x via 

Gi(x)x = gl (x) and G 2 (x)x < g 2 (x) (7) 

where Gi(x) = Gi(x)S and Gi(x) = Gi(x)S. We might also expect the buyer to add additional linear 

constraints, such as requiring the latest shipment from any supplier to arrive earlier than a certain date, or 

requiring all deliveries to arrive the same day. There can also be constraints specific to particular suppliers, 

e.g., the buyer doesn't want any more than 100 units from supplier 5. These can be handled simply as 

constraints on the individual x 1; * and added as extra rows to Gi(x), G 2 (x), gi(x), and i 2 (x)- 

With aggregation, the quadratic form to be minimized is (Ex - mM)^; 1 (x) (Ex - p(x)) subject to the 

constraints given in Eq. (7). This minimization can be carried out through a straightforward generalization 

of the method given in Appendix D. 
i* We can also relate x and x by x = Ex + 1> for some constant vector v % which will not cause any complications. However, 

there seems to be little practical reason to do so. 
16 With no additional computational complexity, we can allow S to depend on x. 
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6.2 Implementation 

In this section we outline an implementation of the entire e-procurement invention. We begin with a high- 
level description of the architecture, then fill in the details by describing a complete object model. 

8.2.1 High-level Architecture of the Invention 

There are at least two modes in which the invention may be used. First, the invention may reside at the 
site of large buyers, and suppliers who wish to sell to the buyer may be required to submit their capabilities 
via a web interface to the buyer. The invention may also be used within a marketplace hosted by a third 
party. Buyers/sellers log onto the market, submit their preference/capabilities, and act on the results. The 
architecture is modular enough to support both modes of operation. 

In Figure 1 we present an architecture for the invention. We describe the architecture, starting from the 
optimization algorithm which finds matches between buyers and sellers and work our way outwards. 

A controller surrounds the optimization engine, feeding it buyer preferences and seller capabilities. If 
multiple optimization processes are running (perhaps on different machines), the controller can also do 
load balancing, forwarding the request to the least busy process. The controller decomposes preferences 
and capabilities into their constituent buyer- and seller-specific versions (see below), selects the most specific 
matching preference/capability pairs, and sends them to the matching engine for optimization. The controller 
then collects responses from the matching engine and returns them to the buyer. Additionally, the controller 
lop all results into a database for recording purposes. 

Another layer, called the Connector in Figure 1, separates the graphical user interface (GUI) through 
which users communicate with the tool from the controller. This layer serves a number of functions. The 
connector transforms the description of preferences and capabilities from the GUI into a form suitable for the 
implementation of the matching engine. Part of this transformation involves validation of appropriate input 
from the GUI layer so that no malformed input is ever sent to the controller. The Connector layer can also 
pull data from ERP or SCM systems and automatically infer preferences (using the total cost of ownership 
function) for the buyer. The enterprise abstraction layer insulates the Connector from the precise details of 
the manner in which the ERP and SCM data needs to be gathered. Total cost of ownership is evaluated in 
the simulation modules, which may either be running locally at the client's site or running centrally at the 
main server. These simulation modules pull operational data (the vector 0)bom the enterprise abstraction 
layer. A preference optimization module (TCO) minimizes the total cost of ownership to determine the ideal 
trade and the flexibilities around the ideal trade. 

At the outmost level, a layer provides integration with the GUI and/or host system. A number of 
administrative systems are expected at this layer. Market administration services allow easy definition 
of trading spaces, the dimensions of negotiation, limits on continuous variables, allowed settings of the 
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discrete variables, etc. User administration services allow an administrator to define buyers, passwords, 
spending limits, etc- Supplier services accomplish similar tasks on the supply side. Managers for preferences, 
capabilities, and match results ensure that these objects are properly stored in a database. This layer layer 
also dynamically generates the html necessary for presentation of the data via a web interface to buyers and 
sellers. 

For maximal portability, communications between the View and Connector are via XML documents. For 
maximal efficiency, communications between the Connector and matching controller are as serialized Java 
objects. 

6.2.2 An Object Model for the Invention 

The fundamental objects required for the invention are preferences from buyers, capabilities from sellers, and 
match results' returned to all parties. The components of such objects have already been considered from a 
mathematical point of view, and we now describe one possible computer representation. 

In this section we describe a complete grammar for the object model. The following syntactic conventions 
are used: 

• (nt) denotes a non-terminal symbol nt 

• [obj] denotes an optional grammar segment obj 

• {obj} denotes 1, or many times the grammar segment obj 

• denotes a production rule for non-terminal symbol If there are multiple rules, say (a), (b), and (c), 
then these are denoted as 

s (nt) -4(a)|(6)|(c). 
In contrast, a production rule of the form 

(nt> -><a>,(&>, (c) 

indicates that the non-terminal (nt) is composed of three grammar segments, (a), (b), and (c) 

• terminal keywords are in serif font 

Obvious non-terminal grammar elements like (string) and (integer) are not described. 
Supply Side 

To represent capabilities that apply to a specific buyer (perhaps for contractual reasons), we have defined a 
capability to be a list of (buyerSpedficCapability). With one exception, a buyer-specific capability applies 
only to one buyer - that buyer associated in the id field of the (buyerSpedficCapability). The exception 
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occura if the id field is * or wildcard. This indicates that the capability applies to all buyers. Using 
buyer-specific capabilities, suppliers can represent specific capabilities to certain buyers and generic capar 
bilities applying to all other buyers. By not including a wildcard (buyerSpecificCapability) and only listing 
(buyerSpecificCapability)s applicable to specific buyers, sellers can also represent the fact that they will 
trade only with a subset of all buyers. In cases where both the wildcard (buyerSpecificCapability) and a 
(buyerSpecificCapability) applicable to a specific buyer apply, the most specific (buyerSperificCapability) is 
selected. 

A schematic of a (sellerSpecificPreference) is given in Figure 2. 
We begin at the top level of a capability: 

capability -4 {(buyerSpecificCapability)} 

where 

(buyerSpecificCapability) -* id: (id), 

discrete: {(discreteVarDescription)}, 
continuous: {(contmuousVarDescription)}, 
range: {(range VarDescription)}, 
[discreteConstraint: (discreteConstraint)], 
instance: {(discreteCapabilitylnstance)} 
[aggregationParticipation: 0|1]. 

(id) identifies a buyer or group of buyers. Individual buyers are represented by some unique identifier. (say 
an integer) and the group of all buyers is identified by the wildcard character **\ So we have 

(id) -4 (integer) | * . 

aggregationParticipation is a Boolean flag giving the supplier's' willingness to participate in aggregate orders 
to the identified buyer. 

Each of the variable constituent components is described by 

(discreteVarDescription) -4 name: (integer), 

allowedValues: {(integer)} 
(continuousVarDescription) -4 name: (integer), 

min: (double), 

max: (double) 
(rangeVarDescription) -4 name: (integer).- 
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la its simplest form, a (discreteConstraint) is a list of more primitive constraints 

(discreteConstraint) {(primitiveDiscreteConstrsunt)} 
where each primitive constraint is composed as follows: 

(primitiveDiscreteConstraint) -> name: (string) 

variables: {(discreteVarName)}, 

includes: 0 1 1, 

values: (integerMatrix) 

(discreteVarName) is the name of the discrete variable involved in the constraint 

(discreteVarName) (integer). 

The includes field is a bit. If the bit is 1, then the combinations listed in the values field are the allowed 
values the variables may take on. If the bit is 0, then the combinations listed in values are the excluded . 
combinations, i.e., everything in the powerset of the variables is allowed except those combinations listed in 
values. The order of the variable names is significant, since they will be assumed to be in the same order in 
values. If there are a variables involved in the constraint, and c constraints, then (integerMatrix) is an a x c 
matrix of integers: 

(integerMatrix) (integer Vector), - • • , (integerVector) 
(integerVector) (integer), • • • , (integer) 

The (discreteCapabilitylnstance) component is described by 

(discreteCapabilitylnstance) -> mask: (discreteVarMask), 

[rangeCapability: {(rangeVarlnstance)}], 
continuousCapability: (continuousCapability) 
continuousConstraints: (continuousConstraints) 

A (rangeVarlnstance) defines a range variable and has the form 

(rangeVarlnstance) -f name: (integer), 
min: (double), 
max: (double). 

The (discreteVarMask) relates to the discussion of 6.2.2. As in Table 3 we have 

(discreteVarMask) -> { (extended Var Value)} 
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where an (extendedVar Value) is either an integer from the domain of the discrete variable or the wildcard 
character '*': 

(extendedVarValue) -» (integer) |*. 

(continuousConstraints) describes the hard linear constraints for the continuous variables. Since these con- 
straints may be either inequality or equality, we have 

(continuousConstraints) -> [equality: (linearConstraints)], 

[Inequality: (imearConstraints)] 

Both the equality and inequality constraints are expressed through a matrix which is c x n c where c is the 
number of constraints, and a vector which is c x 1. Consequently we have 

(linearConstraints) -> matrix: (doubleMatrix), 
vector: (doubleVector) 

A (doubleMatrix) is defined by 

(dqubleMatrix) -» (doubleVector), • • * , (doubleVector) 
and a (doubleVector) is just what the name suggests - a vector of doubles: 

(doubleVector) (double), • • • , (double). 
The only remaining undescribed element above is (continuousCapability) whose description is 

(continuousCapability) breakpoints: (doubleListMatrix), 

valAtBreakPoints: (doubleListMatrix) 

(doubleListMatrix) describes a n c x n c matrix whose elements are lists of (double): 

(doubleListMatrix) -> (doubleListVector),- , (doubleListVector) 
(doubleListVector) (doubleList), — , (doubleList) 
(doubleList) -^{(double)} 

It is assumed that the rows and columns of the matrix are in some canonical order so that we know which 
continuous variable is referenced. A natural order is the one defined in {(contmuousVarDescription)} 

Preferences 

Just as capabilities may be buyer-specific so too may preferences be seller-specific. The same rules deter- 
mining which seller-specific preference to apply are followed. A schematic of a (sellerSpedficPreference) is 
given in Figure 3. 
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We define a preference as follows 

(preference) {(sdlerSpedficPreference)}[, (aggregatedPreference)) 

i.e., a preference is a list of (sdlerSpecificPreference) with an optional aggregated preference. We first 
describe (sellerSpecificPreference) and then consider (aggregatedPreference). 
The (sellerSpecificPreference) is composed as follows 

(sdlerSpetificPreference) ->td: (id), 

discrete: {(discreteVarDescription)}, 
continuous: {(continuousVarDesoription)}, 
range: {(rangeVarDescription)}, 
dimensionWeights: (dimensionWeights) , 
discreteTradeoff: (tradeoffTables) 
[discreteConstraint: (discreteConstraint)], 
instance: {(discretePreferencelnstance)} 

Of these elements, only (dimensionWeights), (tradeofiTables), and (discretePreferencelnstance) have yet to 
be defined. (dimensionWeights) gives the weights of all dimensions that indicate their importance. For 
convenience we break up the weights according to the three types of variables. Thus we have 

(dimensionWeights) -> range: (double Vector), 
discrete: (double Vector), 
continuous: (doubleVector) 

A (doubleVector) has been described previously. Each of the corresponding vectors is as long as the number of 
range, discrete, or continuous dimensions. (tradeofiTables) is an naiscrete * ^discrete matrix of (tradeofFEable): 

{tradeofiTables) -> (tradeoftTableMatrix) 
(tradeofflfcbleMatrbc) -> (tradeofflkbleVector), • • • , (tradeoffTableVector) 
(tradeoflTableVector) -» (tradeofiTabie), ■ • • , (tradeofiTable) 

A (tradeoffEable) is simply a matrix of double values. 

Finally, we turn to the last undefined component of a (preference). A (discretePreferencelnstance) is 
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composed as follows: 

(discretePreferencelnstance) -> mask: (mask), 

[rangeldeal: {{range Varlnstance)], 
continuousldeal: (doubleVector), 
tradeoff Matrix: (doubleMatrix), 
[continuousConstraints: (continuousConstraints)] 

The rangeldeal and continuousldeal fields give the desired range and continuous trade parameters. The trade- 
offMatrix field gives the positive definite matrix expressing the tradeoffs amongst the continuous variables. 
(continuousConstraints) have been described previously in the sell-side specification. 

To complete the specification of preferences, we conclude with the definition of (aggregatedPreference). 
Refer to the discussion of section 6.1.7 for details. 

(aggregatedPreference) -> participants: {(aggSpecification)}, 

contributionType: (contributionTypeVector), 
additionalConstraints: (continuousConstraints), 
discrete: {(discreteVarDescription)}, 
continuous: { (continuous VarDescription)}, 
range: {(range VarDescription)}, 
dimension Weights: (dimension Weights), 
discreteTradeoff: (tradeofETables) 
[discreteConstraint: (discreteConstraint)], 
instance: {(discretePreferencelnstance)} 

In the above definition, the previously defined elements maintain their meaning. The additionalConstraints 
field allows the buyer to express constraints around the aggregation, such as "all orders must arrive on the 
same day," etc. participants lists the suppliers who can participate in the aggregation and their characteristics. 
Note that if the wildcard supplier participates, the order can potentially be aggregated across all suppliers. 
(aggSpecification) describes information specific to a supplier participating in the aggregation. It is defined 
by 

(aggSpecification) -¥ id: (id). 

id identifies the participating supplier and constraints specific to that supplier defined in an accompanying 
(sellerSpecificPreference) will be used in the optimization. Additional information may be added as required. 
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UlOvlCwC lllOOn 


output 


specificity 




{continuousl, rangel} 


0 


jx x * *] 


{continuous2, range2} 


1 


fx! * X 3 ] 


{continuous3, range3} 


2 



Table 3: Example of discrete masks for specifying continuous and range variables which are dependent on 
discrete variables. x x and x 3 signify specific values for the first and third discrete variables. The specificity 
of each mask is indicated in the third column. 

The contributionType field is used to define the f vectors used in aggregation. The {contributionTypeVector) 
consists of n c elements indicating the type of aggregation for each continuous dimension: 

(contributionTypeVector) -> (contributionType), • • • , (contributionType). 

Possible contribution types include 

(contributionType) -4 sum, average, zero. 

sum sets f = 1, average sets £ = l/|£(x)|, and zero sets £ = 0. 
Masking 

We have allowed constraints, ideal values, tradeoffs, and continuous capabUities to be dependent on discrete 
variables. In this section we describe an efficient manner in which to encode this dependence. 

The data structure must represent continuous and range variables for all valid discrete inputs. An 
efficient way to do this is to use hierarchical definitions. At the top of the hierarchy are the definitions 
of the continuous and range variables for the discrete values x* = [*,-••, *]. These values apply to all v 
unless more specialized masks are defined. A more specialized mask of the continuous and range variables is 
specified by defining values for some of the components x*. The more components that are defined, the more 
specialized the definition. The most specific mask is always used. An example definition for three discrete 
variables is given in Table 3. The response to the lookup [xi x 2 £3] is {continuous3, range3} if and only 
if xi = xi A x 3 = xa, {continuous2, range2} if and only if xx = x x A x 3 ^ x 3 , and {continuousl, rangel} 
otherwise. 
Match Results 

Returned to the buyer is a list of matches with different suppliers, which can be ranked and viewed in many 
different ways in the GUI. A (matchResultList) is a list of matchResult: 

(matchResultList) -> {(matchResult)}. 



26 



WO 01/45007 PCT/US00/33017 

A match result may either be a (singleSupplierMatchResult) or an (aggregatedMatchResult): 

(matchResult) -> (angleSupplierMatchResult)| (aggregatedMatchResult). 

A (singleSupplierMatchResult) represents the best trade with a single supplier and is composed of the 
following elements: 

(singleSupplierMatchResult) -» supplierld: (integer), 

utility: (double), 
feasible: 0(1, 
[costFactors: {(double)], 
continuous: {(double)}, 
discrete: {(discreteVarDescription)}, 
range: (range Varlnstance) . 

The supplierld indicates the supplier sourcing this trade and the utility field indicates the utility of the trade 
(which can be used to rank the trades), feasible is a bit indicating whether or not a feasible trade with this 
supplier was found. The continuous, discrete, and range fields list the respective trade parameters determined 
by the matching algorithm. The optional cost factors field lists the constituent costs contributing to the 
total cost of ownership C 0 evaluated at the trade point returned in the (singleSupplierMatchResult). 

An (aggregatedMatchResult) returns the optimal trade when the buyer has requested aggregation. It is 
composed of the following elements: 

(aggregatedMatchResult) -> utility: (double), 

feasible: 0|1, 
[costFactors: {(double)], 

supplierTradePara meters: {(supplierTradeParameters) }. 

As before, the utility field gives the utility of the aggregate trade, and the feasibility flag indicates whether or 
not a feasible aggregate trade was found (there may not be if the constraints were too stringent). costFactors 
can also be returned in C c is sufficiently general to handle aggregated trades. Finally, (supplierTradeParameters) 
lists the trade parameters for each supplier involved in the aggregation. It is defined as follows: 

(supplierTradeParameters) -> supplierld (integer), 

continuous: {(double)}, 

discrete: {(discreteVarDescription)}, 

range: (rangeVarlnstance). 
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6.3 Summary 

We have described an efficient computational procedure in which to encode buyer's trading preferences 
and hard constraints, supplier's delivery capabilities and constraints, and optimize to find those matches 
between one buyer and one or many sellers that maximize the buyer's utility. By optimizing against both 
qualitative and quantitative factors, and exploiting the trading flexibilities possessed by both parties, the 
system determines better trades. The tool is particularly useful as companies move their direct material 
purchasing online. By optimizing across flexibilities, win-win trades are discovered for both trading parties. 

The representation of trading preferences is designed to be expressive yet easily elicitable from a buyer, 
and computationally tractable. The representation of supplier capabilities was chosen to parallel the manner 
in which suppliers already think of their delivery capabilities and seamlessly includes volume discounts and 
incremental costs. These supplier capabilities may be part of an online catalog. The representation of the 
negotiation space is rich, supporting three types of variables. 

We have outlined a manner in which preferences may be inferred automatically through models of the 
purchasing company. Such models incorporate many cost factors, taking the total cost of ownership into 
account. The system provides trades which minimize the total cost and represent significant new savings. 

The invention can operate both at a buyer's site, where suppliers input their capabilities through an 
HTML interface to the world wide web or as an embedded part of an electronic market hosted by a particular 
web site. The invention may operate at regularly scheduled intervals or sporadically in lieu of current request 
for quotations (RFQ). The buyer may broadcast a RFQ event to suppliers, indicating a time within which 
suppliers must respond. At the close of the event, the buyer can use the present invention to assist in the 
analysis of the supplier responses. 

Complex algorithms have been specified which should permit most matching optimization to occur in 
near real-time. The rapidity of optimization, combined with graphical what-if tools, allows for analysis and 
exploration of trades, which should significantly improve the quality of purchasing decisions. 

While the above invention has been described with reference to certain preferred embodiments, the scope 
of the present invention is not limited to these embodiments. One skilled in the art may find variations of 
these preferred embodiments which, nevertheless, fall within the spirit of the present invention, whose scope 
is defined by the claims set forth below. 
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What is claimed in the present invention is 

1. A method for discovery of trades between one or more buyers and one or more sellers comprising the 
steps of: 

(a) expressing one or more terms of an ideal trade and one or more flexibilities by at least one of said 
buyers; 

(b) expressing one or more capabilities by at least one of said sellers; and 

(c) determining at least one optimal trade with respect to said one or more terms and said one or 
more flexibilities of said at least one buyer and said one or more capabilities of said at least one 
seller. 

2. A method for discovery of trades between one or more buyers and one or more sellers as in 1, wherein 
said one or more terms comprise one or more members of the group consisting of continuous factors, 
discrete factors and range factors. 
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A Distance Function for Range Variables 

The distance between buyer range variable n - fcj.n) and seller range variable r, = fo.r,) is related to 
the overlap defined as 

overlapfo.rj) = / dx M{x\ri)M(x\rj) 

where M{x\ r,) is a Gaussian in x centered at m = + with standard deviation a f = a(r i - r { ) and 
a is a tunable parameter. The range of integration is bounded by the seller's capabilities. Evaluating the 
integral, we find 



overlap^,, r,) = , exP 
2^27r(a? + <^) 



> 



By adjusting a we can control the hardness of the boundary at the ends of the range (smaller a corresponds to 
sharper dropoff outside the range). The maximum overlap will occur when ~ m and <tj = 0. Substituting 
those values in the above equation, we get 

maxOverlap = — 

^27r(a? + a?) 

The distance <f(r<,rj) is related to the overlap so that it varies from 0 (for identical range variables) to 
oo (for maximally distant range variables). Consequently, we define 

B Prom Rankings to Numerical Values 

The numerical values in the tables Z(s€ if jcj) can be formed from the buyer's rankings in a number of ways. 
Since we are minimizing distance, we want preferred options to contribute little Z and progressively less 
favorable options to have progressively larger Z values. To approach this, we first consider the problem of 
defining a Z 1 value normalized to lie between 0 and 1 and then take Z = - ln[l - Z% 

One possibility for Z' is a linear scaling. If <k is the integer indicating the number of values . in p i9 
r i°H x i) indicates the buyer-specified rank (smaller rank means more preferable) for the option >q = a € V* 
conditional on x,, then we might set Z'(x* = a 9 x,) = (rj 00 ^) - 1)/(* - 1). Alternatively, we might use 
a non-linear scaling and set Z(x* = a,v t ) = (1 - p#*M- l )f{l - P di ' 1 ) for some 0 € [0, 1] (exponential) 
or set Z(vi = a,vj) = ((rJ o) (wi) - l)/(Jli - 1))" for some n > 0 (algebraic). The corresponding scaling 
functions for Z are given in Table 4. 
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type 

linear 
exponential 
algebraic 



In 



In 



Ri-1 



Ri 



1-0* 



(fi). 

-1 



0r< a >( VJ j-l 



Table 4: Possible numerical ranking functions. 

For those allowed values of x\ which are left unranked by the buyer, we can assign them a common 
numerical value greater than the lowest ranked value. For example if Ti{xj) is the largest ranked value we 
might assign all remaining values to have = n(xj)/(di - 1) in the linear case or = 

(1 _ - or Z(x u *j) = (1 - V(ri(*j) 4- !))/(! - 1/*) in the non-linear cases. 



C Determination of Scaling Factors 

In this section we supply the details necessary to determine the scaling factors Q r and Q d * In the equations 
defining the average distances we have need of the integrals f v du exp[- (u — u^C" 1 (u — u)] and f v du (u - 
aye 1 (u-u) exp[- (ii-u^C 1 (u-u)]. In general, these integrations are not analytically tractable; however 
we will assume that the region of integration defined by x and x extends at least 2 standard deviations out 
in all directions so that the volume of integration is effectively infinite. In this case we can evaluate the 
integrals in closed form: 

J du exp[-(u - uJ'CT^u - u)] « yj- 

and 



7T ne 

detC" 1 



J da (u - u)<C- J (u - &) exph(u - u^C" 1 (u - ft)] « ^ 



7T n * 



detC" 1 * 

To this point we have not been specific about what is meant by the average over the range variables. In 
fact, we will assume that, however we define this average, it is independent of what the buyer requested and 
thus can be removed from inside the sum over the discrete variables. In this case we will concentrate on 
setting the discrete and continuous components equal. We have {d c )/n e — {dd)/n d where 

(dd) _ Qd Ev^exp[-Q d ck] J v du exp[-d c ] ^ (d c ) _ 1 £ v exp[-Q d <k] J v d u d c exp[-d e ] 
n d ~ n d £ v exp[-Q d cfc]/ v du exp[-d c ] n c 

and cancelling denominators Q& is determined by solving 



«c Y, v exp[~Qddd} f v du exp[-d c ] 



^Ytddtxpl-Qddd] f du exp[-d c ] = ±J2*M-Qdd d ] f dud c exp[-d c }. 
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This can be simplified further given our assumptions about the Gaussian integral: 

Qd d d exp[-Q d d d ) = 1 y> exp[-Q d d d ] 
n* V Vdet C-^v) 2 ^ ^detC Hv)' 

In most cases there will also be at least one finite root, which is what we need. If there are a multiplicity 
of roots, we select that root having the smallest value to ensure that the distribution defined by the utility 
function is as flat as possible. If there are no real roots, then we select the Q d which minimizes the error. 

D Phase 1 Optimization 

In this section we detail the solution to the phase 1 optimization problem defined in Eq. (5). 

D.l General Formulation 

We begin by writing down the objective function: 

di(x, x) = (xW -/i^'c-^xW + i!( r;x ) 

where (see Eq. (8)) 

x« = (I - MM(xW, {Or'c^W, {6 W }) {6 |*) } ) + (i _ M «(vW {6iS}))- l M(«({6g})uW (8) 

as long as - x i* > - x **^(*<2 + Tnis must be optimized over x subject to the constraints 

Gi(x)xW = gi(x) and G 2 (*)xW = g2(*) for a given x. Since x is fixed throughout the optimization 
K(r; x) is just a constant which can be evaluated and then ignored as far as determination of x is concerned. 
We simplify our notation so that d x (u) = (xW - pj'C" 1 (xM - p) + iZ(r). Similarly, 

x« = (I-M«({^^ 

where we have defined the n c -vector s = (I - MW)~ l c and the n c x n' c matrix S = (I — MW)" 1 !^ 6 ). 18 
n' c < n e is the number of continuous dimensions upon which the supplier response depends. Putting these 
pieces together we have 

d 1 (xW) = xW t Q({6g,{6g})xW + q«xW (9) 
where Q is an n' c x n' c matrix given by 

Q(fl8.ttB}) - S*({6W,{6iS})C^S({6W,{68}) 

16 In the next section we consider S further, taking advantage of any redundant variables to simplify the QP task. 
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q is an n'-vector given by 

q({*$,{i$» = {©W Wi(J,{aB}) -/«)'• 

and the constant term is 

const = (s({6« {&«}) - -tfcz 1 (8({6iJ, {bg}) -») + Jl(r) 

It is easily shown that the matrix Q is positive semidefinite, d\ (uW) must be optimized subject to the 
constraints 

GiS({6g, {&g})* (t) = & " Gi8({jW, {6g», (10) 
G a S({6«, {6g})x<» < g 2 - G 2 8({6j»>, {&$}), (11) 

^(0<"({*B.{lB}) + S({>{3.{ij!i|})^ ) £^>(0 VI, (12) 
(t) < xW < x<») (t) V t (13) 

where xM(i)' = [ Xl (bg) ... *„(6«)] and x«(i)* = [^(ig + 1) ... + 1)] are the vec- 

tors of breakpoints for the f$ (k e [l,n e ]) functions. Similarly, for the buyer variables, xW(i)* = 

[*i ) *n(&i3)] and = [*i(Ci+l) ••• x„(6|5 + 1)] are the vectors of breakpoints 

for the /JJ (fc e [1, n c ]) functions. 

di is optimized in two steps. First we assume a value for the integer variables {b\'l} and {b^l) fixing the 
constraints in Eqs. (10) and (11). We then minimize the quadratic di subject to these constraints using a 
standard quadratic programming (QP) algorithm 17 to find x^({6$}» {****})• This result is then used to 
find x* («, {&#}, {jjj}) = s({6g}, {$}) + S({6g}, {bg})x^({6|J}..{6g})- By searching exhaustively 
over {b\*l} and {b\*l} we determine 

%«(*) = argmin {*8».*) 

Before turning to the manner in which this optimization is accomplished in section D.3, we discuss a few 
improvements to the input into the QP algorithm. 

D.2 The Quadratic Program 

In this section we discuss improvements in the specification of the quadratic programming task to make it 
better conditioned and more efficient. 

It is possible that the matrix Q used as input to the QP is singular or near singular, depending on the 
rank of Si. Consequently, it is important to use an interior point QP algorithm for which this is not a 
problem. 

"Numerous third parties supply quadratic programming packages, some based on interior-point methods which are appro* 
priate, in the present context where matrix Q may be singular. 
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D.2.1 Variable Elimination 

There may be rows or columns of S which are entirely zero. It is important to take advantage of this structure 
to eliminate unnecessary computation. 

If row i is zero, this indicates that supplier variable u\ a) is independent of any buyer variable and is simply 
equal to If there are any constraints involving only supplier variable t, then these can be checked before 
calling the QP solver to eliminate wasteful computation. 18 If the constraint is satisfied, we can eliminate it 
and then call the QP function. If column j of S is all zero, this means that there is no dependence of u< a > 
on uf* and so there is no need to optimize over that buyer variable. We will, however, still have to check 
that the constraints on are satisfied. This observation allows us to simplify the quadratic programming 
problem by eliminating this unnecessary variable. 

We formulate the problem to take advantage of both these observations. Assume there are m < n c 
non-zero rows of S and m' < n' c non-zero columns. Relabel the rows and columns so that 



Si 0 
0 0 



where Si is of size mxm' and, with a slight abuse of notation, the O's indicate matrices of the appropriate 
size. In like fashion we write 



x<<> 

4* 



£<»> = 



4 ] 



where x£ fi) is an m vector, is an n c - m vector, x[ 6) is an m* vector, and x£ 6) is an n' c - m' vector. 
x[* } indicates those components of xW which vary with buyer variables while x^ indicates those variables 
which are constant. xj 6) are those buyer variables which affect at least one supplier variable and x£ 6) are 
those buyer variables which do not affect any supplier variables. If we also partition s as s = [s{ s'Jthen 

x[ a) = Sl + Six[ 6) and 4 8) = s 2 . 



Finally, we will also have need of the partitioned co variance matrix 



where (C^)' = C^.„ since C" 1 is symmetric. 

Including the dependencies on fi, the quadratic form we must optimize is 



WiCr.lW, W 1 C 1 ^W 2 
WaCj-jWi W^Wa 

18 Additionally, the QP doesnt like to have a row of zeros in the constraint matrix. 



sx-Mx + Siui" 

«2 —1*2 
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where we have made the obvious definitions 



w = 


Wi 


0 


and n = 






0 


w 2 







Expanding the quadratic, the constant term is 

(si - Mi^WiCrj Wi ( Sl - Mi) + (s 2 - M2)*WaC^W 2 (Ba - A ) + 2( Sl - Mi) t W 1 C 1 ^W 2 (s 2 - /x 2 ). 
The linear term is 

q - ((si - MiJ'WiCr.iW^' +2(s 2 -/i^WaCjjWxSOxW 
and the quadratic term is 

Q^^J'SiWiC^WiSiX^. 
The constraints involving the optimization variables uf* are 

Gl;lSiXp* = gi - Gi ; iSi - 6i ;2 8 2 , 

G 2; iSix[ fr) < g 2 - G 2; isi - G 2;2 s 2 , 
where we have simplified notation by dropping the dependence on breakpoints and defined 

G '=h. g S:2 ] 

i.e., G< ; i is the columns of Gi involving the active buyer variables while G< j2 is formed from the inactive 
buyer variables. Similarly, the constraints which can be checked prior to even solving the QP are 

xi a) (0<s 2 <x ( 2 ff) (i) Vi. 

D.2.2 QP Normalization 

There are two types of normalizations we perform for increased numerical stability. If the scales of the 
continuous dimensions vary greatly (i.e., one dimension has an ideal value of 10~ 2 and another dimension 
has ideal value 10 6 ), then the matrix C" 1 becomes badly conditioned, which makes S badly conditioned and 
thus makes the QP unstable. Similarly, if the breakpoints have radically, different lengths, the QP problem 
can also become unstable. We consider how to remedy each of these problems. 

Dimension normalization 
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The quadratic form to be minimized is (x - j*)*cr l (x - /i). We modify this problem by defining the 
dimensionless variables <j>% through = fam in terms of which 

(x - f*)*C~ l (x-fi) = «> - l)*C- l (* » 1) 

where C" 1 = DC~ 1 D and D = diagO*). When making this variable transformation we must assure that 
the constraints on x are adjusted so that Gjx = gi and G 2 x < g 2 become the constraints 9i<j> = gi and 
9i4> < g2 on 0 where £1,2 = Gi (2 D. If there are global bounds of the form x < x < x , then these are 
translated to 0 where £ = D -1 u and = D" 1 !!. 

In addition to these changes, we also need to modify the capabilities appropriately. Recall Eq. (2) which 
reads (after dropping the dependence on the discrete variables) 

k 

In terms of the dimensionless variables we have 

+ (15) 

Thus, when we use the dimensionless variables, the function /,*,* needs to be modified so that the y dimension 
is divided by m and the x dimension is multiplied by ji*. 

Once the optimization has determined an optimal 0, we return to the original variables through x = T)<f>. 

Breakpoint range normalization 

For larger problems it is sometimes important to normalize the input variables for the QP solver. This is 
carried out in order to improve the QP search performance in situations where there is large variation in 
the range of the segments <f${bk) < 4^ < + 1) a&d the magnitude of their values. This can be 

accomplished by introducing the normalization variables = (0< 6) - ^£ 6) (6*))/(4 6) ( b * + l ) " 4 6) 
Here ffi will have minimum value of zero and maximum of one. We can now express <pW as 

*«=A$ W +p (16) 

where A is the diagonal matrix A = diag(^ frl+1) - ~ >^ m+1) - ^ m) ) and p* = 

m V- iw J- 

Expanding Eq.(9) using the definition of £ in Eq.(16), we obtain 

<*i(w; {&*,*}) =|^ (6)t ADQ({6 ilfc })DA# (6) + q'({fc, *})DA*' + p'DQ({6 <(Jb })DA* (6) + 
ip*DQ({6 i(fc »Dp + q^i^BDp. 
Hence the quadratic programming problem to be solved in terms of £ is 

QP(x<>>;{& a }^ (17) 
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where Q = ADQ({6,,*})DA and q' = (q*({6 <(fc })+P t I>Q({& j , fc }))DA. Similarly, the equality and ioequality 
constraints of Eqs. (10) and (11) expressed in terms of ^ are 

6i({M)* (6) = ix *"» 6 2 ({6 i(fc })^ ) <t2, (18) 

where 

G,({6 i , fk }) = G 1 S({b ilfc })DA, 

G 2 ({6 i ,*}) = G 2 S({6 i . fc })DA, 
gi = gi - G lS ({6 llfc }) - G 1 S({6 j .*})Dp, 
ta = g2 - G 2 s({6 ilfc }) - GiS({6 ilfc })D P . 

D.3 Optimization Over The Breakpoint Variables 

We have detailed how to determine x^(x, {b\ % {&$}) given a specification of all {bg} and {&$}, but 
how do we determine the best feasible choice of {&$} and {bg}? Before we answer that question, we have 
to examine the infeasibilities that might arise due to improper choices of «nd 

D.3.1 Elimination of {to,*} Infeasibilities 

Most of the settings of {&$} or will not yield a feasible solution because Eq. (13) or Eq. (12) will be 
violated. The reason is that, for any k, it is highly likely to get two segments x l £\b itk ) < x[ b) < xf> (6 ilfc + 1) 
and x k b) (b jlk ) < x[ b) < a4 k) (6j,* + 1) that do not overlap (i.e., the optimizer will not be able to find a xj^ 
value that satisfies both bounds, and hence no feasible solution will be found). The same is also true for 
supplier response segment x^. 

To eliminate those infeasibilities, we construct a super set of breakpoints {t37|} such that 

*P = IK4^8 = i).»P(*B - so. - .«P<« - «8» 

Hence each piecewise linear function /,-,* will be described by a set of breakpoints xv k = {** U)» • • * » x k (**)} 
and function values /<,ft(x*(l)).--- ,/*,*(**("*))■ These function values at the new breakpoints are lin- 
early interpolated from the old breakpointrfunction pairs. After this reduction we have /t fc - 1 different 
b k € [1, Kk - 1] values, which eliminates the heed for constraint Eq. (13). • 

The continuous supplier variables and constraint Eq. (12) are more problematic. Suppressing any 
dependence on we know 

ft« * 
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If {b\*l} and {b\*} 7 - • • , are known, then this expression allows us to bound the values may take 

on, and thus bound the interval in which the feasible b\*} lie. If we define the functions 



sraax(x,y) if$>0 
sraax(s, x, y) = \ and smin(s, z, y) = < 

3min(x,y) ifs<0 



smin(x,j/) ifs>0 
$ max(x, y) if 5 < 0 



then the minimum and maximum values that x\** may take on are: 

min(*«) = • • ,6{ J_i>. {©) + £«to(«JJ (ig), -PflBUM? + D) 

+ £sm^S(04^©>^ + D). (19) 

and 

max(xW) = c^J,- • ■ ,*g-J, {.*$}) + E^K'i^S).^^),^^,* + 1)) 

+E-«(«8(»B).*f ) (»8).«? , («a-+i)) ( 2 °) 

respectively. The above equations imply that we can fully translate a complete specification of buyer variable 
breakpoint intervals to a corresponding set of feasible supplier variable breakpoint intervals. We will search 
over the breakpoints using a branch-and-bound algorithm, and will need to be able to infer the restriction on 
the values of feasible supplier breakpoint intervals given a partial specification of buyer breakpoint interval 
variables. This inference, which narrows the domains of the supplier breakpoint variables, is called con- 
straint propagation. Before describing how the branch-and-bound algorithm search over breakpoint interval 
variables is accomplished, we show how the constraint propagation may be obtained. 

Narrowing of Supplier Breakpoint Domains Through Propagation 

Let fi)jP be a vector of integers specifying a set of breakpoint intervals for the fcth buyer variable. = 1 if 
the interval is uniquely determined, and \0k\ = k$ - 1 if no restrictions on the interval have been specified. 
We index the /th component of 0^ as 0^ n I € [1, - 1]. A set of breakpoints for all buyer dimensions 
is written as {0^}. In analogous fashion, we define 0%\ 0$, I e [1,4* } " ^ ^ for the supplier 
breakpoint interval variables. 

from a specification {0^} we infer the feasible domains {/^} using max(ar^) and min(x^) where: 

max(xW) =£ max ,{«,*({&<»>}) + smax^M^^WUV + D)} + 

j§^fi->{*' fc({ ^ (21) 



39 



WO 01/45007 



PCT/US00/33017 



and 

min(xW) =£ min + sm^^^ + 1))} + 

£ min UAib?}) + smin(mW(6W),4»)(6W), a: W(6« + 1))} (22) 

We begin with i = 1 (where no supplier variables contribute to x[ 9 ^) and determine the interval I\ = 
[min(x[ ,) ) > max(x ( 1 ' ) )]. We then form the vector with components {j3<^ ^ 0}. 

Equations (21) and (22) can then be iterated to find the interval h — [min(x^ ) ),max(ar^ ) )] using the 
previously calculated min(rc[^) and max(x[**). With J 2 known, the components of 0^ are determined as 
above. Iterating this procedure, we can find the potentially narrowed domains of all {ffi}* 

We now return to the optimization problem: we need to search over all Yik K k? x 11* K $ combinations of 
breakpoints to determine the combination that minimizes d\ (x op t (x, {b^ } , {&j^}), x) . The simplest manner 
in which this accomplished is via exhaustive enumeration over all {b^} and {6^}. However, once all b^ 
have been specified, we need only search exhaustively over the corresponding {&($)} inferred as above. 
There are, however, many ways to speed up the search over the buyer breakpoint intervals. We describe two 
methods. 

D.3.2 Local Search 

This method is extremely fast but not guaranteed to return the global optimum across all buyer breakpoint 
intervals. Here, initial random and feasible 19 {fcj^}> and {bj^} sets are first created. After the continuous 
optimization is carried out using the initial {6^}, and {fcjjf*}, a new {6^}, and {1$} set is generated from 
the old one by examining whether each generated ar^'({6^},{6^}) solution is on the boundaries of the 
interval x^(b k ) < x 1 ^ < x { ^{b k + 1) or x^\b k ) < x^ < x[ 8 \h + 1). If xf ] is on the lower bound of the 
interval, 6£ 6) = b$ - 1; similarly, if it is on the upper bound, b^ = i£ 6) + 1. If xf> is not on either end, 
then - b%\ Of course this done while maintaining Vjp in the interval 1 < &£ fr) < - 1. Supplier 
breakpoint variables are updated in a similar manner. The procedure is repeated until no decrease in distance 
is attained. 

D.3.3 Branch-and-Bound Search 

Branch-and-bound is an efficient procedure used to exhaustively search a space and eliminate configurations 
that can be proven to be suboptimal. The theory behind branch-and-bound is well understood [5]. The 
idea is to structure the space as a tree of partially-completed possibilities and to explore the tree so that 
good complete solutions are constructed early on. Branches of the tree that can be shown to yield solutions 
"Feasible in the sense that they do not violate the constraints in equations (10), (11), (12), and (13). 
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which cannot be better than the currently best solution need not be explored. Heuristics indicate the order 
in which to explore the tree and provide a lower bound on the fitness of partially completed solutions for 
pruning of suboptimal branches.. 

In the present context, we carry out branch-and-bound search over the space of possible buyer breakpoint 
combinations. The tree is structured beginning from a node at which none of the bj^ have been assigned. 
From this node, the tree branches into the - 1 possible assignations for b\ b K From any of these nodes, the 
tree branches into the itf* possible values for The tree branches at each node to assign the remaining 
variables. At the leaves of the tree, all buyer breakpoint variables have been assigned. 

Prom any node, the search proceeds by choosing the most promising node to consider, as indicated by 
a variable ordering heuristic and considering the possible values for that variable in the order indicated 
by the value ordering heuristic. A complete solution and its accompanying distance is obtained when the 
search ends at a leaf of the tree. To avoid unnecessary exploration, a lower bound heuristic indicates when a 
partially completed solution (node) need not be expanded because all completions will have higher distance 
than the best located so far. We describe these three heuristics for search over buyer breakpoints. 

Lower Bound 

Given a set of buyer breakpoint interval vectors {0^} t we can propagate to {jff^} and the intervals 
[min(x^),max(a;^)]. A lower bound on the value di can be obtained by ignoring the supplier capabilities 
and minimizing 

*(x) = (x<'>«/i) t C- 1 (xW-/i) (23) 

subject to 

Gix«= gl (24) 
G 2 x< 8 > < g 2 (25) 
m\n(x\ 8) ) < x\ 9) < maxfr^ ), (26) 

The lower bound thus obtained is denoted by lb({0W}). 
Variable and Value Ordering 

The lower bound lb({/?(*)}) can be used to define an ordering over variables and the values any variable may 
take on. We define the functions L^ t i which operate on the set ifi^} by locking the kth buyer breakpoint 
variable to At any node we can determine the lower bounds lb(Lkj{0W}). We can select to explore 
the variables in order of increasing lb* where 
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i.e., the average over all possible completions for &£\ M Within any given variable the (say k*) values can be 
explored in order of increasing lb(Ljk. ( i({/J (6) })). If we are willing to do more computation, we can explore 
the branches in an interleaved order according to all the lb(L* ( i({/JW})). 

Pruning and Search 

During the branch-and-bound search, branches can be pruned if the node has lb(Lfc,i({0 w })) greater than 
the distance of the lowest examined leaf node. 

Search begins at the root node where each 0^ = [!,*••, - 1], We applying the ordering heuristics 
above to move to a node where one breakpoint variable has been fixed. Irom that node, we move down the 
tree, fixing a value at a time until we reach a leaf. We then progressively backtrack up the tree, eliminating 
branches whose lower bounds are greater than the best leaf examined. 

Once the optimal set of buyer breakpoints has been determined, we search exhaustively over the supplier 
breakpoints feasible with this choice according to Equations (19) and (20). 

Improvements 

We have not yet exploited the fact that, from any node having {j?^}, the constraints in Equation (12) 
may narrow the domains of the buyer breakpoint variables deeper in the tree. The computational cost of 
narrowing the breakpoint domains may be warranted if the pruning thus obtained is sufficient. 

E Phase 2 Optimization 

In this section we detail an efficient algorithm to optimize Eq. (5): 

arg min d!(x(x),x) + <fc(x) subject to C< 6 >(x) AC< a) (x). 
Our approach is a marriage of constraint programming and local search. 

E.l Local Search over ^-neighborhoods 

We begin with a definition. The d-neighbors of a discrete configuration x are all those configurations, x*, 
which differ from x in d attributes. Attribute i differs in two configurations if x* ^ x{. Thus the distance 
varies from 0 if the configurations are identical, up to if all attributes are different. The ^-neighborhood 
of x is the union of all d = 1 neighbors, d = 2 neighbors, up to d = D neighbors. The P-neighborhood 
around x is indicated by Md{*)- We note that, for any configuration, its neighborhood is the entire* 
discrete space. 

Using the D-neighborhood, we apply a standard hiliclimbing algorithm which respects all the discrete 
constraints. D is specified before the algorithm begins and serves as an efficient manner to trade off running 
^We might also take the minimum over all possible completions \h k {{fi^}) = min f <») . lb(L Ai t({0 W }))- 
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Table 5: Configurations examined in the D = 3 neighborhood search around (a) the initial configuration 
11111 and (b) the initial configuration 0+0+0 (duplicate search configurations are indicated in bold) . 

time with the quality of the final solution. The hillclimbing algorithm takes as input a starting configuration 
x t and finds the configuration x t+ i € -A/i>(*t) which minimizes di(x(x),x) + subject to fiW(x) A 
CW(x). The algorithm then iterates again starting from The hillclimbing algorithm terminates when 
it finds a configuration xt which has a distance less than or equal to all of its JD-neighbors. If D = n*, 
this is the global minimum. Given a function which exhaustively searches a D-neighborhood for the lowest 
distance feasible configuration, a J9-neighborhood hillclimbing algorithm is easily implemented. 

There is one noteworthy point about the implementation of an efficient D-neighborhood hillclimber. 
Consider a D = 3 search over n<f = 5 discrete variables. Without loss of generality, we label the initial 
configuration as 11111. In the D = 3 neighborhood search, the configurations examined are given in Table 
5(a) as indicated by the *s. Let us assume that the best neighbor in the D neighborhood of 11111 differs in 
components 2 and 4. We indicate this best neighbor as 1+1+1, which serves as the starting configuration 
for a new D = 3 search. The configurations examined in the second neighborhood search are given in 
Table 5(b). Ftom Table 5(b) we see that the search branches 2, 7, and 10 are redundant, having already 
been performed in the D neighborhood search around lllll. 21 Thus these branches may be pruned in the 
second D neighborhood search. A search branch covers a previous search and may be pruned if the vector 
of positions of the wildcard characters covers the vector of positions at which the previous best solution was 
found. The binary relation covers is defined as 

x covers x* 4=^ x{ € x V xj € x*. 
"The duplication of search effort is most acute when D is large and the fittest neighbor is found very close by (small <Q. 
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The covers relation is transitive but not reflexive. Thus, to increase the efficiency, we can keep a list of the 
previous branches searched and ignored search brandies that have already been covered, 

E.l.l Improvements 

Throughout the present discussion we have assumed that D remains constant through the hillclimbing 
process. There are good reasons why this should not be the case. With any reasonable landscape correlation 
structure, we know that lower distance neighbors are more likely to be found early in the hillclimbing process. 
Thus a reasonable heuristic is to start with a larger D early on and gradually lower it with time. 

E.2 Exhaustive Search of a ^-neighborhood 

To complete the discussion of the phase 2 optimization, we describe a computationally efficient way in which 
to find the lowest distance feasible configuration within a D-neighborhood. We use a branch-and-bound 
method with heuristics tuned for the present problem. The factors tuned include the variable order in which 
we search, the order of the values any given variable may take, and a heuristic which says when certain 
branches of the search space may be safely ignored. 

We use heuristics which early on in the search generate low-distance configurations, and later prune parts 
of the search that we can prove must have higher distance than the best configuration generated so far. 

E.2.1 Heuristics And Lower Bounds 

We focus on obtaining good heuristics for the discrete contribution to the total distance, then later incorpo- 
rate the continuous and range contributions since these are simpler. 

At any point in the search, some variables have been assigned and some have not. We divide up the sum 
over the Zij{xu Xj) into four contributions 

i locked f locked < locked j floating 

i floating j locked i floating j floating 

Since we are minimizing and all Z itj (x u *>i> and Qd are positive, we can prune any search node which 
has Z\ t \ > than the distance of the lowest distance configuration found thus far, since additional terms will 
only add to the sum. 

We obtain a lower bound on the distance by lower-bounding each of the above four terms. Just taking 
the lowest value in the tables can yield differing values for any given variable, e.g., the values for a single 
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variable x*: 

x;=argmin ^(xs.x,), x? =argmin Zj.^xi), and x? = arg min Z M (xi) 
may all be different. To correct this we write 

t floating 

where and Zf,f are calculated as before and 

j locked 

The best choice for the floating variable x< is 

xf-argminZftW (27) 

In this way all the possible choices for x* are at least consistent, and the true lower limit is still guaranteed 
to be higher than this (since the Z f(f contribution can only be higher). 

This approach also naturally suggests a variable ordering and a value ordering within any variable. We 
can try the variables in increasing order of their corresponding (*f) 80 * hat wW>le iy which has the 
smallest Z$(xJ)> is tried first. For any given variable x* we can also order the possible value assignations 
for x< according to the array of values so that the value x* == xj is tried first and the value which 

maximizes Z$(xi) is tried last. 

With little extra work we can impose a consistency requirement raising the lower bound for Z f(f . We 
decompose the sum as follows 

Zf.f = Qd 52 WiW&A****) = Qd 52 VHWj{Zij{xt*Xi) + 

floating:/^ t,i floatiiig:i<i 

Instead of taking the minimum component-wise as previously suggested, t.e instead of 

Zf t t -Qd 52 j? in [ w i w 3 Z *J (**> » 
< l Jfloating:Jyi ,^<r,^, 

we take the minimum so that 

Zti = Qd 52 minf^^Cx^x^ + Z^Cx^Xi))]. (28) 

ij floating: j<f 

which enforces pairwise consistency between variables i and j. 
The lower bound, lb, we have obtained is thus 
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where is determined from Eq. (27) and Z£ f is determined from Eq. (28). 
Integration Including Continuous and Range Contributions 

In this section we summarize the complete heuristic specification, including the effects of the continuous and 
range variables. As above, we use the definitions 

Z\ f \ = Qd 52 w * w i^u(*i» x j). 

ij locked 
j locked 

^f*f = Qd 53 min[ttf t n; i (Zf (i (x <1 x 5 ) + Z;,i(**, x*))] 
IJ floatlng:j<* 

so that a lower bound to the complete discrete contribution to the distance is Z\ f \ + Zfo 4- £ f 2£f(*T) where 
^rrargniin^^^xi). 

We now extend this result to include the continuous and range variables. The continuous and range 
variables depend on the discrete variables through the mask. Since we expect there will be relatively few 
masks we exhaustively enumerate all the discrete possibilities. Consequently, we define di(xfl 0 ating|*locked) 
as the distance contribution made by the continuous and range variables (the d\ contribution to Eq. (5)). 
We define the continuous + range contribution for any given floating variable i as 

dW(Xi) S min di(Xfloatlng blocked) 
flo&tlng\i 

where the min is taken over all floating variables other than i and where floating variable * is locked to value 

Since each continuous/range/discrete distance contributes additively to the total distance, we can rank 
the possible values for any given floating variable i according to Z$(xi) + <4 f) (i/<). The best (according to 
this heuristic) setting for floating variable i is 

xf =argmln{^J(^)+c4%)}. 

The heuristic then suggests a variable ordering in order of increasing Z$(xf) + ^(vj) so that we first 
select the floating variable with the lowest sum and last select the variable with the highest sum. 

Finally, we turn to the lower bound on the distance. Noting that the continuous and range distance 
contribution does not have the form of a sum over the dimensions (as the discrete contribution does), we 
write the lower bound, lb, as 22 

where i* = arg mini (^ff (*?) + <^( v i)) 0- e -» ** labeb the dimension with the least 2$ + 4° 

contribut ion to the distance). 

is the Kronecker delta function which is 1 if and only if * = at*. 



46 



WO 01/45007 PCT/USOO/33017 
E. 2 .2 Improvements 

We point out that an extensive literature exists on parallel branch and bound algorithms. Many of the 
techniques developed are applicable to the present problem. 
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AMENDED CLAIMS 
[received by the International Bureau on 8 May 2001 (08.05.01); 
new claims 3-108 added; remaining claims unchanged (14 pages)] 

1. A method for discovery of trades between one or more buyers and one or more sellers 
comprising the steps of: 

5 (a) expressing one or more terms of an ideal trade and one or more flexibilities by at 

least one of said buyers; 

(b) expressing one or more capabilities by at least one of said sellers; and 

(c) determining at least one optimal trade with respect to said one or more terras 
and said one or more flexibilities of said at least one buyer and said one or more 

10 capabilities of said at least one seller. 

2. A method for discovery of trades between one or more buyers and one or more sellers 
as in 1, wherein said one or more terms comprise one or more members of the group 
consisting of continuous factors, discrete factors and range factors. 

3. (New) A system for determining one or more trades between a buyer and one or more 
15 suppliers comprising: 

(a) one or more variables defining a space of negotiation; 

(b) a utility function of said one or more variables for expressing a utility of the one or 
more trades to the buyer over the space of negotiation comprising: 

i. an ideal trade to the buyer defined by one or more ideal values corresponding 
20 to said one or more variables; and 

ii. at least one flexibility in at least one of said variables expressing how the utility 
of the trade to the buyer varies in the space of negotiation for said ideal trade; 

(c) one or more capabilities for defining a subspace of the negotiations space wherein 
the one or more suppliers have an ability to trade; and 

25 (d) an optimizer determining at least one of the trades that is optimal with respect to 

said utility of the buyer subject to the capabilities of the one or more suppliers. 

4. (New) A system for detennining one or more trades between a buyer and one or more 
suppliers as in claim 3 wherein said one or more variables comprise one or more of the 
following: 

30 (a) one or more continuous variables x, one or more discrete variables h and one or 

more range variables r. 

5. (New) A system for detennining one or more trades between a buyer and one or more 
suppliers as in claim 4 wherein each %i of said one or more continuous variables x has an 
allowed range over which said each continuous variable a* may vary Xi 6 Xi = 

35 wherein £ is a lower bound of said continuous variable Xi and Xf is an upper bound of 

said continuous variable Xi\ 
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6. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 5 wherein each x» of said one or more discrete variables x has a 
value from a domain 

Hi € Vi = [1, ■ • ■ »*1 where <k > 0 is an integer giving the number of possible values 
5 that said discrete variable x$ may assume. 

7. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 6 wherein the space of negotiation comprises a tensor product 
Xi ® • • ■ ® X ne ® 2?i ® • • • ® T> nd . 

wherein 

10 n c is the number of said continuous variables; 

na is the number of said discrete variables; 

Xi is said allowed range of said continuous variable Xf } and 

Di is said domain of said discrete variable x$. 

8. (New) A system for determining one or more trades between a buyer and one or more 
15 suppliers as in claim 4 wherein said utility function u((x, x, r)) comprises an expression 

of a distance function d(x, x, r) that defines a distance from said ideal trade in the space 
of negotiation. 

9. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 8 wherein said distance function comprises at least one of the 

20 following: a continuous distance, a discrete distance Z(x), and a range distance JL 

10. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 9 wherein said range distance depends on the value of at least one 
of said discrete variables x. 

11. (New) A system for determining one or more trades between a buyer and one or more 
25 suppliers as in claim 10 wherein said range distance is defined as: 

«(r;x) = ££ x iii(ri;*) 
wherein: 

r is an nr vector of tables of preferred values of said one or more range variables, = 
30 n r is the number of said range variables; and 

n > u- 

12. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 10 wherein said range distance is a distance dfarj) between said 
range variable r< = (U^i) of the bu y er 8111(1 said ran 6 e r i ~ (Hi> f i) o{ at least 

35 one of the suppliers. 
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20 



13. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 12 wherein said distance between said buyer range variables and 
said supplier range variable d(r f> r,) comprises an overlap between said buyer range 
variable and said supplier range variables, overlap (r<, r,). 

14. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 13, wherein said overlap is defined as 



overlap^, r,-) = / dx M{x\ u)M (x; r 3 ) 



where 

10 M{x\rj) is a Gaussian in x centered at (ij = (tj + *j)/2 with standard deviation Oj = 

orfo - Lj)\ 

M(x;ri) is a Gaussian in x centered at m = (u+^i)/ 2 with standard deviation <n = 

ot(ri-u); 811(1 

a is a tunable parameter. 

15 15. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 14 wherein said range distance is defined as 

vi \ i„ r °verlap(ri,rj) 1 
i?(r i ,r 3 ) = -ln^ max0veriap J 



wherein 



erf[jg] 

maxOverlap = 



16. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 9 wherein said continuous distance is quadratic and is determined by 
a positive semi definite n c x n c matrix C" 1 wherein n c is the number of said continuous 
variables. 

25 17. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 16 wherein said continuous distance is defined as 

(x-/i(x)) t C- 1 (x)(x-/iW)- 

wherein 

30 fi is an n c -vector of ideal values of said continuous variables that may depend on said 

discrete variable x. 
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18. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 9 wherein said continuous distance is a convex function. 

19. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 9 wherein each of said discrete variables *i has a value from a 

5 domain D{. 

20. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 19 wherein said discrete distance Z(x) maps a discrete space V x ® 
. . . ® V nd onto the positive real line [0,oo] wherein n d is the number of said discrete 
variables. 

10 21. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 20 wherein said discrete distance Z{x) is defined as 



2(*)=E{2i(*i)+ E 3u(**.**)} 
i=i 



wherein 

15 each Zij is a table comprising djdj entries; and 

didj is the distance if said ith discrete variable has value x it conditional on said jth 
discrete variable having value Xj. 

22. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 21 wherein values in said tables Z { j are determined from one or 

20 more rankings by the buyer. 

23. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 22 wherein said discrete distance Z is defined as Z = - ln[l - Z% 

wherein 

Z 1 is normalized to lie between 0 and L 

25 24. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 23 wherein said normalized distance Z l is a linear scaling. 

25. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 23 wherein said normalized distance Z' is an exponential scaling. 

26. (New) A system for determining one or more trades between a buyer and one or more 
30 suppliers as in claim 23 wherein said normalized distance Z' is ah algebraic scaling. 

27. (New) A system for determining one or more trades between a buyer and one or more 
supplies as in claim 9 wherein said discrete distance Z{x) is a function of one or more 
pairs of said discrete variables X{ } Xj. 
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28. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 8 wherein said distance is generated from a ranking of preferred 
values for said one or more variables by the buyer. 

29. (New) A system for determining one or more trades between a buyer and one or more 
5 suppliers as in claim 4 wherein said utility function u((x,x,r)) expresses one or more 

tradeoffs among the one or more variables x, r. 

30. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 8 further comprising one or more scaling factors, Q c , Qd> Qr to nor- 
malize contributions of said at least one continuous variable x, said at least one discrete 

10 variable x, and said at least one range variable r to said utility function u((x,x*,r)) for 

defining a baseline wherein the one or more variables contribute equally to said utility. 

31. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 30 wherein said distance function with said normalized contribution 
is defined as: 

15 d(x, x, r)) = d c + Q d d d + Q T dr 

wherein 

d c is a contribution to said distance by said continuous variables, 
Qddd is a contribution to said distance by said discrete variables, and 
Q r dr is a contribution to said distance by said range variables. 

20 32. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 30 wherein values of said scaling factors are set so that average 
distances of said one or more from variables from said corresponding one or more ideal 
values are equal. 

33. (New) A system for determining one or more trades between a buyer and one or more 
25 suppliers as in claim 32 wherein said average distance comprises utility-weighted average 

distances over said space of negotiation for placing more weight on better ones of the 
trades. 

34. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 33 wherein said utility- weighted average distances are defined as 

E„ J v du Q r d r txp[-Qrd r - Q d d d - d e j _ £ M expf-Q^dal S, ttxpt~Q r d r )/y <fu cxp[-d e ] 
6%} * r ' ~ E r Sy du **P[-Qrd r - Qddd ~ dc) E„ expt-Qdd,,] exp[~Q^d r ) J v du exp{-d e J 

= Em Er Sy dn Q d d d cxp[-Q r d,. - Q d d d - d e ) _ d d npi-QM^ exp[-Q r dr-) J v du exp[-d c j 

Em E r Sy exp[-<? r d r - Q d d d - dc} E* **p[-Qdd d } £ r exp[-Q r d P ] J v du exp[-d e ] 

= Zm E r Sv d» d c exp[-q r d r - Qddd ~ d e ) = exp[-Q a d rf ) £ p exp{-Q r dr.) J v du d c exp[-d c ) 
EvEr/v^ «cp[-Q r d r - Q d d d - d c ) JZ M npi-Qddd) En «p(-Q r d r } J v du exp[-d c J 

wherein 

35 Ylx indicates a repeated sum £ Xl • • * £ over ^ possible discrete trades; 

£) r indicates a sum over all of said range variables; and 
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Q c = 1 because Q r is interpreted as Q r /Qc and Q d is interpreted as Qd/Qc- 

35. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 34 wherein said scaling factors Q Cy Q r are determined from 
said utility-weighted average distances. 

5 36. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 30 further comprising one or more weights to enable the buyer to 
weight said contributions of said one or more variables to said utility. 

37. (New), A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 36 wherein said distance function comprises at least one of the 

10 following: a weighted continuous distance, a weighted discrete distance Z w (x) and a 

weighted range distance. 

38. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 37 wherein said weighted continuous distance is defined as: 

15 (x-MMj'ctfMfx-^x)). 

wherein 

C-^WcC^Wc; 

W c is a diagonal matrix formed from w c ; 
w c is an n c -vector of weight for said continuous variables; and 
20 fi is an n c - vector of ideal values for said continuous variables that may depend on said 

discrete variables x. 

39. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 37 wherein said weighted discrete distance Z w (x) is defined as: 

25 wherein 

Zij(ni, Xj) is the distance if said ith discrete variable has value >q conditioned on said 
jth discrete variable having value xj; and 

u/dji is the ith component of a n<* vector of weights for said discrete variables. 

40. (New) A system for determining one or more trades between a buyer and one or more 
30 suppliers as in claim 37 wherein said weighted range distance is defined as 

wherein rir is the number of said range variables, 

r is an fir-vector of tuples of preferred values of said one or more range variables, 
n = and 
35 n > 
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41. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 37 wherein said weighted distance function comprises at least one of 
the following: an n c -vector of weights for said continuous variables, w c , an n^-vector of 
weights for said discrete variables, and an iv-vector of weights for said range variables. 

5 42. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 41 wherein said weights are normalized so that 

c- 1 = W C C^W C 

43. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 41 wherein values for said weights depend on values of said discrete 

10 variables, >q. 

44. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 4 farther comprising a total cost of ownership function expressing 
a total cost of membership to the buyer that varies over the negotiation space. 

45. (New) A system for determining one or more trades between a buyer and one or more 
15 suppliers as in claim 44 wherein said total cost of ownership function comprises one or 

more cost contribution, 

46. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 45 wherein said one or more cost contributions comprise one or 
more of the following: piece part costs, freight costs, setup costs, quality assurance 

20 costs, repair costs, and revenue generated from the trade. 

47. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 44 wherein said total cost of ownership function is defined as 
C o (x,x,r;0) 

wherein 

25 f3 represents one or more other factors comprising one or more of the following: forecasted 

demand and current inventory levels. 

48. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 47 wherein said one or more other factors are extracted from at 
least one of the following: an enterprise resource planning system and a supply chain 

30 management system. 

49. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 48 wherein said one or more other factors are extracted in real time 
for enabling continuous, real time optimization. 

50. (New) A system for determining one or more trades between a buyer and one or more sup- 
35 pliers as in claim 47 wherein minimization of said cost of ownership function C 0 (x, x, r; /J) 

determines said ideal trade to the buyer: Xo P t(/J), *opt08), r op t(0) 
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51. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 50 wherein said at least one flexibility is determined by a Hessian 
matrix H = [hij] 
wherein hij is defined as 



dxidxj 



x=x o pt(0),x=xopt(0),r=roi>t(0) 



52. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 3 wherein said one or more capabilities comprise one or more of the 
following: price discounts on large volume orders, and variation in delivery as a function 

10 of price. 

53. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 8 wherein said one or more capabilities specify one or more of the 
following: 

one or more continuous capabilities, one or more discrete capabilities, and allowed values 
15 for said range variables. 

54. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 53 wherein said allowed values for said ranges variables contribute 
to said distance function. 

55. (New) A system for determining one or more trades between a buyer and one or more 
20 suppliers as in claim 53 wherein said allowed values for said range variables comprise one 

or more pairs (Ej,fj) wherein Vj is a lower bound for the jth one of said range variable 
and Tj is an upper bound for the jth one of said range variable. 

56. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 52 wherein said continuous capabilities are one or more responses 

25 . from said suppliers to a request for the buyer. 

57. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 56 further comprising a vector- valued function f(x^ 6 \x( 5 \x) to 
determine said one or more supplier responses wherein 

x^ is said buyer request and 

30 x( 5 ) is at least a previous one of said supplier responses. 

58. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 57 wherein said vector- valued function f(x^,x( s ),x) comprises one 
or more components fa for corresponding ones of said continuous variables: 

*W-/,(x<»,xM) 
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59. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 52 wherein said one or more components, are piecewise linear 
functions. 

60. (New) A system for determining one or more trades between a buyer and one or more 
5 suppliers as in claim 59 wherein said one or more components, /» are specified with one 

or more breakpoints. 

61. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 53 wherein said one or more discrete capabilities of said suppli- 
ers comprise one or more constraints from said suppliers on said one or more discrete 

10 variables. 

62. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 3 wherein said one or more capabilities are represented by one or 
more piecewise linear functions. 

63. (New) A system for determining one or more trades between a buyer and one or more 
15 suppliers as in claim 4 further comprising one or more constraints involving said one or 

more variables which must be satisfied for the buyer and at least one of the sellers to 
trade. 

64. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 63 wherein said one or more constraints comprise one or more of 

20 the following: discrete constraints for expressing one or more allowed and/or disallowed 

combinations of values for said discrete variables x and continuous constraints for setting 
one or more requirements on said continuous variable x. 

65. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 64 wherein said continuous constraints from the buyer are linear. 

25 66. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 64 wherein said continuous constraints comprise at least one of 
inequality constraints and equality constraints. 

67. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 64 wherein said continuous constraints depend on values of said 

30 discrete variables, x. 

68. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 63 wherein said one or more constraints comprise one or more of 
the following: 

required delivery time, and an unacceptable color. 

35 69. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 64 wherein at least one of said continuous constraints depend on 
values of at least one of said discrete variables x. 



56 

AMENDED SHEET (ARTICLE 19) 



WO 01/45007 



PCT/USOO/33017 



determining values of said continuous variables that minimize said distance for one or 
more settings of said discrete variables. 

79. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 78 further comprising the steps of: 

5 representing said distance by a function of said discrete variables; and 

determining an optimal one of said settings of said discrete variables by minimizing said 
function of said discrete variables. 

80. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 78 wherein said determining values for said continuous variable step 

10 is performed under one or more constraints on said continuous variables. 

81. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 79 wherein said representing said determining an optimal one of 
said settings of said discrete variables step is performed under one or more constraints 
on said discrete variables. 

15 82. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 4 further comprising means for aggregating at least one of the 
suppliers to participate in said trade with the buyer. 

83. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 82 wherein said aggregating means perform the steps of: 

20 .determining one or more subsets of said suppliers that satisfy one or more constraints 

on said discrete variables. 

84. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 83 wherein said one or more discrete variable constraints comprise at 
least one of the following: buyer discrete variable constraints and seller discrete variable 

25 constraints. 

85. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 83 further comprising the step of optimizing over said continuous 
variables to determine an optimal one of said subset of buyers. 

86. (New) A grammar for a system that determines one or more trades between a buyer and 
30 one or more suppliers comprising: 

(a) one or more capability rules representing one or more capabilities of said suppliers 
to trade with the buyer; 

(b) one or more preference rules representing one or more preferences of the buyer; and 

(c) one or more match rules representing matches between said one or more capabilities 
35 and said one or more preferences. 
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70. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 3 wherein said determining at least one of the trades that is optimal 
step comprises the steps of: 

selecting one of said suppliers; 

5 determining at least one of the trades that corresponds to a maximum value of said 

utility of the buyer and that is within the capabilities of said selected supplier. 

71. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 70 wherein said determining at least one of the trades that is optimal 
step further comprises the steps of: 

10 selecting another of said suppliers; and 

repeating said determining at least one of the trades step for said another selected 
supplier. 

72. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 71 wherein said determining at least one of the trades that is optimal 

15 step further comprises the steps of: 

choosing at least of the suppliers having the highest said maximum value of said utility 
of the buyer. 

73. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 72 further comprising a subsystem to perform said determined trade 

20 between the buyer and the chosen supplier. 

74. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 3 further comprising a subsystem to perform said determined trade. 

75. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 3 wherein said utility comprises at least one of the following: quan- 

25 titative factors and qualitative factors. 

76. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 4 wherein said determining at least one of the trades that is optimal 
step further comprises the step of: 

minimizing a distance from said ideal trade. 

30 77. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 76 wherein said distance comprises one or more of the following: a 
continuous component, a discrete component, and a range component. 

78. (New) A system for determining one or more trades between a buyer and one or more 
suppliers as in claim 77 wherein said minimizing a distance from said ideal trade step 
35 comprises the steps of: 
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87. (New) A grammar for a system that determines one or more trades as in claim 86 
wherein said one or more capabilities of said suppliers are specific to particular ones of 
said buyers. 

88. (New) A grammar for a system that determines one or more trades as in claim 86 wherein 
5 said capability rules comprise one or more of the following: 

(a) a discrete variable rule for representing a description of a discrete variable; 

(b) a continuous variable rule for representing a description of a continuous variable; 
and 

(c) a range variable rule for representing a description of a range variable. 

10 89. (New) A grammar for a system that determines one or more trades as in claim 88 wherein 
said description of the continuous variable comprises at least one of a minimum value 
and a maximum value for the continuous variable. 

90. (New) A grammar for a system that determines one or more trades as in claim 88 
wherein said one or more capability rules comprise one or more constraint rules repre- 

15 senting constraints on value of at least one of said discrete variables, and said continuous 

variables. 

91. (New) A grammar for a system that determines one or more trades as in claim 89 
wherein said one or more constraint rules comprise at least one matrix for representing 
said constraints. 

20 92. (New) A grammar for a system that determines one or more trades as in claim 90 wherein 
said constraints on said values of said discrete variables comprise one or more permitted 
value continuations for the discrete variables. 

93. (New) A grammar for a system that determines one or more trades as in claim 88 
wherein said range variable description comprises at least one of a minimum value and 

25 a maximum value for the range variable. 

94. (New) A grammar for a system that determines one or more trades as in claim 90 wherein 
said constraints on said values of said continuous variables comprise one or more of the 
following: an inequality, an equality, a linear constraint, and a non-linear constraint. 

95. (New) A grammar for a system that determines one or more trades as in claim 86 
30 wherein said one or more capability rules further comprise an aggregation flag indicating 

a willingness of the supplier to participate in an aggregation for the buyer. 

96. (New) A grammar for a system that determines one or more trades as in claim 86 wherein 
said preferences of the buyer are specific to at least one of said suppliers. 

97. (New) A grammar for a system that determines one or more trades as in claim 86 wherein 
35 said preference rules comprise one or more of the following: 
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(a) a continuous variable rule for representing a description of a continuous variable; 

(b) a discrete variable rule for representing a description of a discrete variable, and 

(c) a range variable rule for representing a description of a range variable. 

98. (New) A grammar for a system that determines one or more trades as in claim 97 
5 wherein said preference rules further comprise one or more weights for representing an 

importance of at least one of said discrete variable, said continuous variable and said 
range variable. 

99. (New) A grammar for a system that determines one or more trades as in claim 97 wherein 
said preference rules comprise at least one of the following: 

10 (a) a first field representing an ideal value for said range variable; 

(b) a second field representing an ideal value for said continuous variable; and 

(c) a matrix representing one or more tradeoffs of said continuous variables. 

100. (New) A grammar for a system that determines one or more trades as in claim 97 wherein 
said preference rules comprise a matrix representing one or more tradeoffs of said discrete 

15 variables. 

101. (New) A grammar for a system that determines one or more trades as in claim 97 further 
comprising at least one aggregation rule comprising at least one of the following: 

(a) a list of one or more of said suppliers that can participate in the one or more trades 
with the buyer; 

20 (b) one or more contribution type fields for specifying contribution types of said or 

more continuous variables; and 

(c) one or more constraints around the aggregation. 

102. (New) A grammar for a system that determines one or more trades as in claim 101 
wherein said contribution types comprise at least of the following: sum, average and 

25 zero. 

103. (New) A grammar for a system that determines one or more trades as in claim 101 
wherein said constraints around the aggregation comprise requiring that all orders arrive 
on the same day. 

104. (New) A grammar for a system that determines one or more trades as in claim 99 wherein 
30 said one or more preferences rules further comprise: 

(a) at least one mask for allowing at least one of said ideal value for said range variable, 
said continuous variable, and said one or more tradeoffs of said continuous variables 
to be dependent on values of said discrete variables. 

105. (New) A grammar for a system that determines one or more trades as in claim 88 wherein 
35 said one or more match rules comprise at least one of the following: 
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(a) a single supplier match rule describing at least one optimal one of said one or more 
trades with a single one of the suppliers; and 

(b) an aggregate supplier match rule describing at least one optimal one of said one or 
more trades with an aggregation of said suppliers; 

5 106. (New) A grammar for a system that determines one or more trades as in claim 105 
wherein said single supplier match rule comprises at least one of the following: 

(a) an identifier for indicating said supplier of said trade; 

(b) a utility for indicating a utility of said trade; 

(c) a feasibility flag for indicating whether a feasible one of the trades with said single 
10 supplier was found; 

(d) a continuous variable field indicating a value for said continuous variable; 

(e) a discrete variable field indicating a value for said discrete variable; 

(f ) a range variable field indicating a value for said range variable; and 

(g) a cost factors field indicating constituent costs contributing to a total cost of own- 
15 ership at said trade. 

107. (New) A grammar for a system that determines one or more trades as in claim 104 
wherein said aggregate supplier match rule comprises at least one of the following: 

(a) a utility field indicating a utility of said trade; 

(b) a feasibility field indicating whether a feasible one of said trades with the aggrega- 
20 tion of suppliers was found; 

(c) a cost factors field indicating constituent costs contributing to a total cost of own- 
ership at said trade; and 

(d) a list of one or more trade parameters for said suppliers in the aggregation. 

108. (New) A grammar for a system as in claim 107 wherein said list of trade parameters 
25 comprise at least one of the following: 

(a) an identifier for identifying one of said suppliers in the aggregation; 

(b) a continuous variable field indicating a value for said continuous variable; 

(c) a discrete variable field indicating a value for said discrete variable; 

(d) a range variable field indicating a value for said range variable. 
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AMENDMENT TIMBER ARTICLE 19 



In response to the International Search Report, please amend the above-identified 
application by substituting the attached sheets 28-41 containing the original claims 1-2 and 
new claims 3- 108 for the previous version of sheet 28 containing original claims 1-2. 

Because substituting attached sheets 28-41 alters the page numbers of the sheets that 
follow the claims (i.e. the list of references, Abstract and appendices), Applicant also attaches 
sheets 42-60, corresponding to the previous version of sheets 29-47. No substantive changes 
were made to the previous version of sheets 29-47. Instead, their page numbers were changed 
from 29-47 to 42-60 respectively to accommodate the inclusion of the substitute claim sheets. 
Accordingly, please further amend the above-identified application by substituting the 
attached sheets 42-60 for the previous version of sheets 29-47. 
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Figure 3: Schematic of a supplier-specific preference with examples indicating potential input. 
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